2201.13296v3 [astro-ph.SR] 25 May 2022 


e 
° 


arXiv 


DRAFT VERSION May 26, 2022 
Typeset using LATEX twocolumn style in AASTeX62 


An Isolated Stellar-Mass Black Hole Detected Through Astrometric Microlensing* 


Kariasu C. Sanu,!}t Jay ANDERSON,’ STEFANO CASERTANO,! Howarp E. Bonp,”)? ANDRZEJ UDALSKI,*: + 
Martin Dominik,*'? ANNALISA CALAMIDA,! ANDREA BELLINI,! THOMAS M. BROWN,’ MARINA REJKUBA,’ VARUN BaJaJ,} 
No& Kains,°' 1 Henry C. Fercuson,! Curis L. FRYER,’ AND PuiLip Yock® 


PRZEMEK MR06z,? SZYMON KozLOwSKI,? PAWEL PIETRUKOWICZ,? RADEK POLESKI,® JAN SKOWRON,”® IGOR SOSZYNSKI,® 
Micuat K. SzyMANskI,® KRzyszTOF ULACZYK,*? AND LUKASZ WYRZYKOWSK?° 
(OGLE CoLLaBoRATION) 


RICHARD Barry,!” Davip P. BENNETT,!”!! Ian A. Bonp,’? Yur Hirao,” STELA Isurrant Siiva,!"!4 Jona Konpo,! 
Naoki Kosuimoto,'!? CLEMENT RANC,! NicHoLas J. RATTENBURY,!° TAKAHIRO SuMI,!? DaISUKE SuzuKI,”? 
Pau. J. TRISTRAM,!’ AND AIKATERINI VANDOROU!: |! 

(MOA COLLABORATION) 


JEAN-PHILIPPE BEAULIEU,'* ! JEAN-BAPTISTE MARQUETTE,” ANDREW COLE,'® PascaL Fouqus,”! Kym Hitt,'® 
STEFAN Dieters,'® CHRISTIAN CouTuRES,’? DiJANA DOMINIS-PRESTER,”2 CLARA BENNETT,”? ETIENNE BACHELET,”* 
JOHN MENzIES,”” MICHAEL ALBROW,”° AND KAREN POLLARD”? 

(PLANET COLLABORATION) 


ANDREW GOULD,?”?° JENNIFER C. YEE,?? WILLIAM ALLEN,’ LEONARDO A. ALMEIDA,?!3? GRANT CHRISTIE,?? 
? ? b ? ? 
JoHN DRuMMonpD,*“"®? AvisHAy GAL-YAM,°° EVGENY GORBIKOV,®” FRANCISCO JABLONSKI,?> CHUNG-UK LEE,”® 
Dan Maoz,*° ILAN MANULIS,*! JENNIE McCormick,’? Tim Natuscu,??*? RIcHARD W. PoGGE,”® AND 
Yoss! SHVARTZVALD“* 
(uFUN COLLABORATION) 
UFFE G. JORGENSEN,’ Kuatip A. ALsuBAI,“° MICHAEL I. ANDERSEN,*” VALERIO BozzA, 
SEBASTIANO CALCHI Novati,’? MARTIN BurGcporr, Tobias C. Hinsp,°*”*? MARKUS HUNDERTMARK,!’ 4 
? ¥ ? ? 
Tim-OLIVER Husser,** EAMONN KERINS,”” PENELOPE LONGA-PENA,’® Luict MANcINI,?"?"°S MarTHEW PENNY,”? 
SOHRAB RAHVAR,°” DAVIDE Ricct,°! SEDIGHE SAJADIAN,” JESPER SKOTTFELT,°? COLIN SNopDGRASS,°* 4 
‘3 + J ? ] ? 
JOHN SOUTHWORTH,°? JEREMY TREGLOAN-REED,°©° JOACHIM WAMBSGANSS,!” AND OLIVIER WERTZ°" 
’ ? ? 
(MINDSTEP ConsorTIUM) 


YIANNIS TSAPRAS,!° RACHEL A. STREET,”* D. M. BraAmicu,°*® Kerra Horne,” ** AND IAIN A. STEELE 
(ROBONET COLLABORATION) 


3 


48,49 


70 


1 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 
2 Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 
3 Astronomical Observatory, University of Warsaw, Al. Ujazdowskie 4, 00-478 Warszawa, Poland 


4University of St Andrews, Centre for Exoplanet Science, SUPA School of Physics & Astronomy, North Haugh, St Andrews, KY16 98S, 
UK 


5 European Southern Observatory, Karl-Schwarzshild-Strape 2, 85748 Garching bei Miinchen, Germany 
6 Department of Physics & Astronomy, Barnard College, Columbia University, 3009 Broadway, New York, NY 10027, USA 
7 Center for Theoretical Astrophysics, Los Alamos National Laboratory, Los Alamos, NM 87545, USA 
8 Department of Physics, University of Auckland, Auckland, New Zealand 
° Department of Physics, University of Warwick, Gibbet Hill Road, Coventry, CV4 7AL, UK 
10 Code 667, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
11 Department of Astronomy, University of Maryland, College Park, MD 20742, USA 
12 Institute of Natural and Mathematical Sciences, Massey University, Auckland 0745, New Zealand 
13 Department of Earth and Space Science, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan 
14 Department of Physics, The Catholic University of America, Washington, DC 20064, USA 
15 Zentrum fiir Astronomie der Universitit Heidelberg, Astronomisches Rechen-Institut, Ménchhofstr. 12-14, 69120 Heidelberg, Germany 
16 Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand 


Corresponding author: Kailash C. Sahu 
ksahuQstsci.edu 


* This research is based in part on observations made with 
the NASA/ESA Hubble Space Telescope, obtained from the Space 
Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS 5-26555. 


SAHU ET AL. 


17 University of Canterbury Mt. John Observatory, P.O. Box 56, Lake Tekapo 8770, New Zealand 
18 School of Natural Sciences, University of Tasmania, Private Bag 37 Hobart, Tasmania 7001 Australia 


19 Sorbonne Universite, UPMC Univ Paris 6 et CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, 
France 


20 Laboratoire d’astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geoffroy SaintHilaire, 33615 Pessac, France 
21 Université de Toulouse, UPS-OMP, IRAP, Toulouse, France 
22 Okayama Astrophysical Observatory, National Astronomical Observatory of Japan, Asakuchi, Okayama 719-0232, Japan 
23 Department of Physics, Massachussets Institute of Technology, Cambridge, MA 02139, USA 
24 Las Cumbres Observatory Global Telescope Network, 6740 Cortona Drive, Suite 102, Goleta, CA 93117, USA 
25 South African Astronomical Observatory, PO Box 9, Observatory 7935, South Africa 
26 University of Canterbury, Department of Physics & Astronomy, Private Bag 4800, Christchurch 8020, New Zealand 
27 Maz-Planck-Institute for Astronomy, Kénigstuhl 17, 69117 Heidelberg, Germany 
28 Department of Astronomy, Ohio State University, 140 W. 18th Ave., Columbus, OH 43210, USA 
29 Center for Astrophysics | Harvard & Smithsonian, 60 Garden St., Cambridge, MA 02138, USA 
30 Vintage Lane Observatory, Blenheim, New Zealand 
31 Fscola de Ciéncias e Tecnologia, Universidade Federal do Rio Grande do Norte, Natal - RN, 59072-970, Brazil 
32 Programa de Pés-Graduagéo em Fisica, Universidade do Estado do Rio Grande do Norte, Mossoré - RN, 59610-210, Brazil 
33 Auckland Observatory, Auckland, New Zealand 
34 Possum Observatory, Patutahi, New Zealand 
35 Centre for Astrophysics, University of Southern Queensland, Toowoomba, Queensland 4350, Australia 
36 Benoziyo Center for Astrophysics, Weizmann Institute of Science, 76100 Rehovot, Israel 
37 School of Physics and Astronomy, Raymond and Beverley Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel Aviv 69978, Israel 
38 Instituto Nacional de Pesquisas Espaciais, Astrophysics Division, Sao Jose dos Campos, Brazil 
39 Korea Astronomy and Space Science Institute, Daejon 34055, Republic of Korea 
40 School of Physics and Astronomy, Tel-Aviv University, Tel-Aviv 6997801, Israel 
41 Department of Particle Physics and Astrophysics, Weizmann Institute of Science, 76100 Rehovot, Israel 
Pp y. phy 
42 Farm Cove Observatory, Centre for Backyard Astrophysics, Pakuranga, Auckland, New Zealand 
43 Institute for Radio Astronomy and Space Research (IRASR), AUT University, Auckland, New Zealand 
44 Department of Particle Physics and Astrophysics, Weizmann Institute of Science, Rehovot 76100, Israel 
45 Centre for ExoLife Sciences, Niels Bohr Institute, University of Copenhagen, Oster Voldgade 5, 1850 Copenhagen, Denmark 
46 Qatar Environment and Energy Research Institute (QEERI), HBKU, Qatar Foundation, Doha, Qatar 
47 Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen, Denmark 
48 Dipartimento di Fisica “E.R. Caianiello,” Universita di Salerno, Via Giovanni Paolo I 132, 84084 Fisciano, Italy 
49 Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Napoli, Italy 
50TPAC, Mail Code 100-22, Caltech, 1200 E. California Blud., Pasadena, CA 91125, USA 
51 Lianenweg 7a, 22529 Hamburg, Germany 


52 Institute of Astronomy, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Torun, ul. Grudziadzka 5, 
87-100 Torun, Poland 


53 Chungnam National University, Department of Astronomy, Space Science and Geology, Daejeon, South Korea 
54 Institut fur Astrophysik, Georg-August-Universitdt Géttingen, Friedrich-Hund-Platz 1, 37077 Gottingen, Germany 
55 Jodrell Bank Centre for Astrophysics, Alan Turing Building, University of Manchester, Manchester, M13 9PL, UK 
56 Centro de Astronomia, Universidad de Antofagasta, Avenida Angamos 601, Antofagasta 1270300, Chile 
57 Department of Physics, University of Rome “Tor Vergata,” Via della Ricerca Scientifica 1, 00133 Roma, Italy 
58INAF — Astrophysical Observatory of Turin, Via Osservatorio 20, 10025 Pino Torinese, Italy 
59 Louisiana State University, 261-B Nicholson Hall, Tower Dr., Baton Rouge, LA 70803-4001, USA 
60 Department of Physics, Sharif University of Technology, PO Box 11155-9161, Tehran, Iran 
61INAF — Padova Astronomical Observatory, Vicolo dell’Osservatorio 5, 35122 Padova, Italy 
62 Department of Physics, Isfahan University of Technology, Isfahan 84156-83111, Iran 
63 Centre for Electronic Imaging, Department of Physical Sciences, The Open University, Milton Keynes, MK7 6AA, UK 
64 Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh, EH9 3HJ, UK 
65 Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK 
nstituto de Investigacién en Astronomia y Ciencias Planetarias, Universidad de Atacama, Copayapu , Copiapo, Atacama, ale 
°6 Instituto de I igacid A fa y Ciencias Pl ias, Universidad de A Ci 485, Copiapé, A Chil 
‘pace Sciences, Technologies, and Astrophysics Researc nstitute, University of Liége, Liége, Belgium 
879. Sct Technologi dA hysics R h (STAR) Insti Uni } Liége, Liége, Belgi 


88 Center for Astro, Particle, and Planetary Physics, New York University Abu Dhabi, P.O. Box 129188, Saadiyat Island, Abu Dhabi, 
UAE 


69 Division of Engineering, New York University Abu Dhabi, P.O. Box 129188, Saadiyat Island, Abu Dhabi, UAE 
70 Astrophysics Research Institute, Liverpool John Moores University, Liverpool, CH41 1LD, UK 


AN ISOLATED STELLAR-MaAss BLACK HOLE 


(Received January 28, 2022; Revised May 19, 2022; Accepted May 24, 2022) 
Submitted to ApJ 


ABSTRACT 


We report the first unambiguous detection and mass measurement of an isolated stellar-mass black 
hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source 
star of the long-duration (tg ~ 270 days), high-magnification microlensing event MOA-2011-BLG- 
191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of 
the Galactic bulge. HST imaging, conducted at eight epochs over an interval of six years, reveals 
a clear relativistic astrometric deflection of the background star’s apparent position. Ground-based 
photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of the Earth’s 
motion on the microlensing light curve. Combining the HST astrometry with the ground-based light 
curve and the derived parallax, we obtain a lens mass of 7.1+1.3 Mo and a distance of 1.58+0.18 kpc. 
We show that the lens emits no detectable light, which, along with having a mass higher than is possible 
for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper 
motion for the BH. The proper motion is offset from the mean motion of Galactic-disk stars at similar 
distances by an amount corresponding to a transverse space velocity of ~45 km s~!, suggesting that the 
BH received a “natal kick” from its supernova explosion. Previous mass determinations for stellar-mass 
BHs have come from radial-velocity measurements of Galactic X-ray binaries, and from gravitational 
radiation emitted by merging BHs in binary systems in external galaxies. Our mass measurement is 








the first for an isolated stellar-mass BH using any technique. 


1. MEASURING THE MASSES OF BLACK HOLES 
1.1. Black Holes in Binary Systems 


Stars with initial masses greater than ~20 Mo are ex- 
pected to end their lives as black holes (BHs) (e.g., Fryer 
& Kalogera 2001; Woosley et al. 2002; Heger et al. 2003; 
Spera et al. 2015; Sukhbold et al. 2016). Objects of these 
masses constitute roughly 0.1% of all stars, leading to 
the expectation that the Galaxy should now contain of 
the order of +108 BHs (Shapiro & Teukolsky 1983; van 
den Heuvel 1992; Brown & Bethe 1994; Samland 1998). 

However, the actual detection of stellar-mass BHs 
is observationally challenging, and determining their 
masses even more so. BHs have been identified in the 
Galaxy and Local Group through X-ray emission due 
to accretion in short-period binary systems, most of 
them soft X-ray transients. In such cases, dynamical 
masses of the BHs can be measured or estimated through 
radial-velocity measurements and light-curve modeling 
for the optical companion stars (the techniques are re- 
viewed by Remillard & McClintock 2006 and Casares 
& Jonker 2014). Masses of nearly two dozen BHs in 
X-ray binary systems have been determined using these 
methods, with varying degrees of precision. These “elec- 


+ PLANET Collaboration, MiNDSTEp Consortium 
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tromagnetically measured” BH masses show a distribu- 
tion peaking near 7-8 Mo, with few if any below ~5 Mo 
(e.g., Ozel et al. 2010; Farr et al. 2011; Kreidberg et al. 
2012; Corral-Santana et al. 2016). This suggests that 
a “mass gap” exists between the lowest-mass BHs, and 
the highest measured masses of neutron stars (NSs) in 
binary radio pulsars of ~2.1—2.3 Mo (Linares et al. 2018; 
Cromartie et al. 2020). Recently, however, a few non- 
accreting or weakly accreting BHs have been discovered 
in longer-period spectroscopic binaries in the field (e.g., 
Thompson et al. 2019; Jayasinghe et al. 2021) and in 
globular clusters (Giesers et al. 2019), lying in the NS- 
BH gap with dynamical masses of ~3-4.5 Mo. Precision 
astrometry of nearby stars by Gaia shows the promise 
of detecting additional wide binary systems containing 
quiescent BHs (e.g., Chawla et al. 2021; Janssens et 
al. 2021, and references therein) and measuring their 
masses. At the high-mass end, the BH mass distribution 
falls off above ~10 Mo, and very few electromagnetic 
BH masses are known above ~15 Mo, the only exception 
in the Milky Way being an updated mass determination 
of 21.1+2.2 Mo for the BH in Cygnus X-1 (Miller-Jones 
et al. 2021). Among extragalactic X-ray binaries, BH 
masses as high as 15.65+ 1.45 and 17+4 Mo have been 
reported for M33 X-7 (Orosz et al. 2007) and NGC 300 
X-1 (Binder et al. 2021), respectively. An even higher 
mass of at least 23.1 Me was reported for the compact 
object in IC 10 X-1 (Silverman & Filippenko 2008), but 
this has been questioned (Laycock et al. 2015). 


‘ 
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The first detections of gravitational waves (GWs) by 
the Laser Interferometry Gravitational Wave Observa- 
tory (LIGO; Abbott et al. 2016) revealed a population of 
massive merging binary BHs, BH-NS pairs, and binary 
NSs at extragalactic distances. In the source catalogs 
from the third observing run of the Advanced LIGO and 
Advanced Virgo collaboration (Abbott et al. 2021a,b), 
the inferred masses of the BHs among the pre-merger 
systems range from 6 to 95 Mo, with two low-mass 
outliers among the secondary components at ~2.6 and 
2.8 Mo, which could be either BHs or NSs. 


1.2. Isolated Black Holes 


The electromagnetic and GW mass measurements de- 
scribed above are all for BHs in binary systems, includ- 
ing those undergoing mass accretion or mergers. How- 
ever, there are reasons to believe that a substantial frac- 
tion of stellar-mass BHs are single, rather than belong- 
ing to binaries. First, about 30% of massive stars are 
born single (see Sana et al. 2012; de Mink et al. 2014). 
Second, in a close binary system, the pair may enter 
into a common envelope and merge before the super- 
nova (SN) explosion (e.g., Fryer et al. 1999; Zhang & 
Fryer 2001; Tutukov et al. 2011; Dominik et al. 2012). 
Lastly, in a wide binary, the “natal kick” imparted to 
the companion by the SN event may be large enough to 
detach the two components, producing an isolated BH 
(e.g., Tauris & van den Heuvel 2006; Belczynski et al. 
2016). Being less altered by interactions with compan- 
ions, single BHs potentially provide a more direct probe 
of BH formation than those in binaries. 

Isolated BHs are extremely difficult to detect directly. 
They emit no light of their own, and the accretion rate 
from the interstellar medium is generally likely to be too 
low to produce detectable X-ray or radio emission (see, 
however, Agol & Kamionkowski 2002; Fender et al. 2013; 
Tsuna & Kawanaka 2019; Scarcella et al. 2021 for the 
case of isolated BHs in dense environments). In fact, 
until now, no isolated stellar-mass BH has ever been 
unambiguously found within our Galaxy or elsewhere. 

Microlensing is the only available method for measur- 
ing the masses of isolated BHs. Astrometric microlens- 
ing—the relativistic deflection of the apparent position 
of a background star when a compact object passes in 
front of it—provides a direct method for measuring the 
masses of BH lenses. High spatial-resolution interfero- 
metric observations of microlensing events (Dong et al. 
2019; Zang et al. 2020), and observations of rare events 
where the lens passes over the surface of the source (Yoo 
et al. 2004), can also yield the masses of BH lenses. In 
this paper, we describe how the technique of astromet- 
ric microlensing is used to determine masses. We dis- 


cuss our ongoing program of astrometric measurements 
of microlensing events with the Hubble Space Telescope 
(HST). Then we report the first detection of an isolated 
BH and our measurement of its mass. 


2. MEASURING THE MASSES OF ISOLATED 
BLACK HOLES WITH ASTROMETRIC 
MICROLENSING 


2.1. Microlensing Events and Black-Hole Candidates 


A microlensing event occurs when a star or compact 
object (the lens) passes almost exactly in front of a back- 
ground star (the source). As predicted by general rel- 
ativity (Einstein 1936), the lens magnifies the image of 
the source, producing an apparent amplification of its 
brightness. The lens also slightly shifts the apparent 
position of the source (Miyamoto & Yoshii 1995; Hag et 
al. 1995; Walker 1995)—an analog of the deflection of 
stellar images during the 1919 solar eclipse (Dyson et 
al. 1920), which provided support for the general theory 
of relativity. 

Microlensing survey programs, including OGLE 
(Udalski et al. 2015), MOA (Bond et al. 2001), and 
KMTNet (Kim et al. 2016), carry out photometric 
monitoring of rich stellar fields in the Galactic bulge. 
These surveys typically detect >2000 events toward the 
Galactic bulge annually. To date, more than 30,000 mi- 
crolensing events have been discovered and monitored 
by these survey programs. 

The characteristic scale of gravitational microlensing 
is provided by the angular Einstein radius 9g, given by 


_ 4G Miens TLS 
‘B= V2 TAU? (1) 


where Miens is the mass of the lens, and 





Ts = ™ — Ts = (1 AU) (s-- a) (2) 


is the relative lens-source parallax, with Dy and Dg be- 
ing the distances from the observer to the lens and to 
the source, respectively. 

The Einstein radius 0g, however, cannot be obtained 
directly from the magnification light curve, whose only 
characteristic that carries a physical dimension is the 
timescale. Given the relative proper motion pips be- 
tween lens and source, we can straightforwardly de- 
fine tj, the time for the source to traverse an angu- 
lar distance of 0g in the barycentric reference frame, as 
th = 6z/uLs (more details in §2.4). The distribution of 
the timescale t%, of the observed events peaks around 25 
days, with t§ ranging from a fraction of a day to several 
hundred days (e.g., Wyrzykowski et al. 2015). 
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If BHs constitute a small but non-negligible fraction 
of the total stellar mass of the Galaxy, as described 
above, then a few of the observed microlensing events 
are expected to be due to BHs. Equation 1 shows that 
the angular Einstein radius 0g is proportional to the 
square root of the lens mass. So, all else being equal, 
events due to massive compact objects would preferen- 
tially tend to be characterized by longer event durations 
(tg 2 150 days) combined with an apparent lack of light 
contribution from the lens. However, a degeneracy be- 
tween lens mass and proper motion remains. Thus a 
long-duration event with no light contribution from the 
lens could arise from a high-mass, non-luminous BH lens 
with a large Einstein radius—but it could alternatively 
be due simply to an unusually slow-moving, faint, low- 
mass ordinary star. If the lens were a luminous massive 
star, this would generally be recognizable through the 
contribution of its light, particularly when observations 
are available in two different bandpasses. 

Indeed, several long-duration OGLE and MOA mi- 
crolensing events have been suggested as being due to 
BHs (e.g., Bennett et al. 2002; Mao et al. 2002; Minniti 
et al. 2015; Wyrzykowski & Mandel 2020). However, 
these claims remain statistical in nature, being derived 
from assumptions about the transverse-velocity distri- 
butions. 

The degeneracy between mass and relative velocity 
can be lifted if precise astrometry is added to the pho- 
tometry of the microlensing event. The size of the 
expected astrometric shift is small—of the order of 
milliarcseconds—but it is proportional to the angular 
Einstein radius 0g, and therefore if this small shift can 
be measured, then the mass of the lens can be deter- 
mined unambiguously, as described in detail below. 


2.2. Photometric Microlensing 


Photometric microlensing is the apparent transient 
brightening that results as a background source passes 
almost directly behind a foreground lens (see reviews by 
Paczyniski 1996, Gaudi 2012, and Tsapras 2018). 

With 6g and 0, denoting the angular positions of the 
source and lens as seen by the observer, we can define a 
dimensionless source-lens separation 


Os — OL 





As the lens intervenes near the line of sight from the 
observer to the source, the gravitational bending of light 
leads to a time-varying magnification 

u2+2 


MOS hee 


which depends solely on u = |u|. This expression holds 
as long as the finite angular size of the source star can 
be neglected, which we will assume in the following dis- 
cussion. We will demonstrate in §9.3 that this is a valid 
approximation for the case analyzed in this paper. 

If Fy is the intrinsic source flux, and Fp the back- 
ground flux contributed by any other objects not re- 
solved from the observed source star, the observed flux 
of the target for a specific telescope and filter is given 
by 


Alu()| +9 


F(t) = FsAl[u(t)] + Fg = Fhase 
1l+g 


» (5) 
where Fhase = Fs + Fp is the baseline flux and g = 
Fg/Fs is the specific blend ratio. 


2.3. Astrometric Microlensing 


Microlensing also produces an astrometric shift of the 
apparent position of the source. If we assume that we 
can observe the centroid of light formed by the images of 
the source without any contribution from other bodies 
such as the lens or other neighboring stars, its shift is 
described by the vector 


5(u) = 


U 


+a 6) 


(Gould 1992; Paczyriski 1998). In contrast to the pho- 
tometric microlensing signature, the astrometric signa- 
ture is explicitly proportional to the angular Einstein 
radius Op. 

Moreover, while the magnification diverges (for a 
point-like source) as u — 0, the astrometric shift be- 
comes maximal for u = V2. The light magnification 
falls rapidly with increasing u: for large separations, 
u > 1, the brightness enhancement, A(w) — 1, falls as 
1/u*. On the other hand, the centroid shift, given by 
Equation (6), decreases more slowly with u, and for large 
separations it falls only as 1/u. The astrometric pertur- 
bation thus has a considerably longer duration than the 
photometric signal. For more details, see Dominik & 
Sahu (2000), Sahu et al. (2014), and Bramich (2018). 

The photometric and astrometric signatures of a mi- 
crolensing event are connected because they arise from 
the same source-lens trajectory, u(t). Specifically, if a fit 
to the light curve of a microlensing event already yields 
u as a function of time, the astrometric data then pro- 
vide a direct measurement of the angular Einstein radius 
Op, as well as the orientation angle of the trajectory. 


2.4. Parallarx Effect and Proper Motion 


For constructing the source-lens trajectory u(t), we 
have to consider the proper motions pg and pry of the 
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source and lens objects, as well as their parallaxes, 7g 
and 7. 

Let y(t) (1 AU) denote the projection of the Earth’s 
orbit onto a plane perpendicular to the direction toward 
the source star. The apparent geocentric positions of 
source and lens star are then given by (cf. An et al. 
2002; Gould 2004) 


Os(t) =Os,0 + (t — to) ts — ts Y(t) , 
1 (t) =OL.0 + (t — to) ML — TL Y(t), (7) 


so that for 6(t) = Og(t) — OL (t), one finds 
6(t) = (5 — O1,)o — (t—to) wus + ™rs V(t), — (8) 


where pps = My — Ms and mis = 7, — Tz are the relative 
proper motion and relative parallax between lens and 
source, while (8g — 01,)o = 93,0 — OL. 

Consequently, with the microlensing parallax param- 
eter Te = Trs/Op, u(t) = O(t)/Om takes the form (cf. 
Dominik et al. 2019) 


u(t) = uo + (t — to) to + TE OY(E), (9) 
where 
ty = (to) = OS 9H)0 5 eey(to), (10) 
E 
tho = ti(to) =- T= + me ¥(to), (11) 
as well as 
y(t) = y(t) — ¥(to) — (t — to) ¥(to) . (12) 


By construction, dy(to) = 0 and d¥(to) = 0. 

By choosing to so that up L to, we can write u(t) in 
its components toward northern and eastern directions 
as 











t—t 
Un(t)= 5 — cosy — to sin + ™ dya(t) 
E 
t—t 
ue(t) = —— sin + up cos + mp dye(t), (13) 
E 
where ug = |uo|, te = 1/|uo|, and w denotes the di- 


rection angle of % measured from north toward east. 
Alternatively, the trajectory can be parameterized as 











t—t 
Un (t) = = cos w* — ug sind” + TE Y(t), 
js aaah tage . t 14 
Ue(t) = in sin ~* + ug cosy* + Tp 7e(t), (14) 





where up = |As(t§) — OL (t)|/On, th = On/|eis|, and ~* 
denotes the direction angle of —pys = ps — PL mea- 
sured from north toward east. Note that the starred 





Figure 1. Source-lens trajectory u(t) as seen by the ob- 
server, showing the effect of annual parallax. At epoch to, the 
tangent to the source-lens trajectory to is not (anti-)parallel 
to the direction of the lens-source proper motion fy, but 
differs by 7p (to) [Equation (11)], related to the orbital ve- 
locity of the Earth at to. Consequently, we distinguish the 
direction angles ~, w*, and yrs, referring to to, and +fy5, 
respectively. Furthermore, tg = 1/|to| and t% = 05 /|m1s|- 





quantities, ~*, tf, tj and uj, refer to parameters in the 
barycentric reference frame, whereas the corresponding 
unstarred quantities refer to parameters as seen by an 
observer on Earth. 

This applies to any orthonormal reference frame; the 
only difference is in the specific angle, e.g., Weq, Weel; 
Yea (and Weg, Vea, Vgai) for equatorial, ecliptic, and 
galactic coordinates, respectively, which are related by 
a rotation of the coordinate axes at the target position. 

While most discussions of photometric microlensing 
events choose an ecliptic coordinate frame, the observed 
astrometric data are more easily described in an equato- 
rial coordinate frame, and therefore we adopt the latter 
in the following analysis. In this frame, 


vis = Veq + 180° (15) 


gives the position angle of the proper motion of the lens 
with respect to the source ps, measured from equatorial 
north toward east. We illustrate the geometry of the 
source-lens trajectory u(t) in Figure 1. 


2.5. Measuring the Lens Mass 


As stated above, the only useful physical parameter 
in a typical microlensing event is the timescale th = 
Oe/t~xLs, which is the time it takes the source to traverse 
the radius of the Einstein ring, which itself depends on 
the lens mass, Miens, and the lens-source parallax, 7s. 

However, for long-duration events, the annual parallax 
tends to lead to prominent departures in the photomet- 
ric signature (Gould 1992; Alcock et al. 1995), so that 
a microlensing parallax parameter mg = 7is/Op can be 
inferred. Coincidently, the BH-mass lenses tend to im- 
ply such long-duration events. On the other hand, the 
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astrometric signature is proportional to the angular Ein- 
stein radius 6g, so that by combining photometric and 
astrometric observations, Miens, ™Ls, and pigs become 
fully decoupled. Specifically, with 7 from the photom- 
etry and 6g from the astrometry, the definition of Op, 
Equation (1), immediately gives us 


On et AU) On 
Miens = = , 16 
: Teh 4G KITE 9) 





where & = 4G'/[c?(1 AU)] ~ 8.144mas M5". 

In the case of microlensing toward the Galactic bulge, 
the source often lies at the distance of the bulge it- 
self, which can be verified from its baseline position 
in a color-magnitude diagram (CMD). If spectro- 
scopic observations are available in addition to baseline 
photometry—as is the case for the event discussed in this 
paper—a more accurate source distance can be deter- 
mined. The lens-source relative parallax, p95 = 7E Op, 
can then be used to estimate the distance to the lens, 
using Equation (2). As a bonus, the event timescale 
gives a direct measure of the relative transverse velocity 
of lens and source (which, for stellar remnants, might 
include “kicks” received in SN explosions). This method 
thus provides independent measurements of three sepa- 
rate physical parameters of the lens: its mass, distance, 
and transverse velocity. 





2.6. Characteristics of Astrometric Deflections 


Some features of astrometric deflections under various 
scenarios are described by Dominik & Sahu (2000). To 
illustrate a typical case, we show in Figure 2 the cal- 
culated astrometric shifts and light magnification for a 
nominal event of a BH lens of mass 5 Mo, at a distance 
of 2 kpc from the Sun, passing in front of a background 
source situated in the Galactic bulge at a distance of 
8 kpc. The closest angular approach is assumed to be 
at a separation of 0.05 6g. In this case, the size of the an- 
gular Einstein ring is 0g ~ 4 mas, so that the maximum 
astrometric shift is ~1.4 mas (occurring at a separation 
of u = V2), and the maximum light magnification is a 
factor of ~20 (at closest angular approach). 

As Figure 2 illustrates, and as discussed in §2.3, the 
duration of the astrometric deflection is considerably 
longer than that of the photometric magnification. This 
makes it necessary to carry out the astrometric mea- 
surements over a longer time interval than the photom- 
etry. Although the deflection measured at any given 
epoch provides in principle an estimate of Op, it is nec- 
essary to observe at multiple epochs in order to separate 
the shifts caused by microlensing from those caused by 
the proper motion of the source; observations at a late 
epoch are particularly useful for this purpose. The fig- 


ure also shows that the astrometric shift is close to zero 
at the time of highest magnification; therefore observa- 
tions near the photometric peak are also very useful to 
constrain the source proper motion. 

In Figure 3, we plot the maximum astrometric shifts 
for a source in the Galactic bulge at 8 kpc, as func- 
tions of lens mass. The lenses are assumed to be located 
at distances of 2 and 4 kpc (“disk” lenses), and 6 kpc 
(“bulge” lenses). The dotted line at the bottom shows 
the nominal astrometric precision of 0.2 mas achiev- 
able with high-SNR HST imaging, as discussed below. 
Therefore the deflection is detectable at lo per epoch 
for lens masses down to ~ 0.5Mo, except at lens dis- 
tances larger than 6 kpc. The most favorable situation 
for a precise mass measurement, of course, would be for 
a nearby, high-mass lens. 


2.7. High-Precision Astrometry 


Although an unambiguous determination of lens mass 
is possible from a combination of photometry and as- 
trometry, as we have just discussed, the expected astro- 
metric shifts are extremely small, of the order of mil- 
liarcseconds or less. HST has demonstrated its capabil- 
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Figure 2. Astrometric shift (top panel) and light magnifi- 
cation (bottom panel) for a microlensing event produced by 
a 5 Mo black hole at a distance of 2 kpc passing in front of 
a background star at 8 kpc. The assumed minimum impact 
parameter is uo = 0.05. The maximum astrometric shift of 
the source is ~1.4 mas, at u = V2, and the maximum magni- 
fication is ~20, at the time of closest angular approach. Note 
the much longer duration for the astrometric shift, compared 
to that of the light magnification. 
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Figure 3. Maximum astrometric shift of a source in the 
Galactic bulge at 8 kpc as a function of lens mass, for lens 
distances of 2, 4, and 6 kpc. Since the maximum astromet- 
ric shift occurs at a lens-source separation of u = V2, while 
the maximum light magnification occurs at the minimum 
value of u, all high-magnification microlensing events will 
pass through the point of maximum astrometric deflection. 
High-S/N imaging with HST allows measurements of astro- 
metric shifts with a precision of ~0.2 mas per observation 
epoch, shown by the dotted line at the bottom. 


ity to carry out sub-milliarcsecond astrometry through 
a variety of techniques. For example, high-S/N HST 
observations of isolated sources were used to achieve 
sub-milliarcsecond accuracy, leading to the measure- 
ment of proper motions for several distant hyperveloc- 
ity stars (Brown et al. 2015). A collection of sources 
was used as probes to achieve an astrometric accu- 
racy of ~12 microarcseconds, in order to measure the 
transverse velocity of M31 (Sohn et al. 2012). Spatial- 
scan techniques have been used to achieve an astro- 
metric accuracy of ~30 microarcseconds (Casertano et 
al. 2016; Riess et al. 2018) in the trigonometric par- 
allax of Cepheids used for accurate determination of 
Ho, and to measure the distance to the globular cluster 
NGC 6397 (Brown et al. 2018). Recently, our group used 
the astrometric-microlensing technique to measure the 
mass of the nearby white dwarf Stein 2051 B, achieving 
an astrometric precision of ~0.2 mas per epoch (Sahu 
et al. 2017). Kains et al. (2017) looked for astromet- 
ric deflections in HST observations of 10 microlensing 
events with timescales of <50 days. They achieved an 
astrometric precision of 0.2 mas per epoch (but did 


not detect any deflections). From the ground, Zurlo 
et al. (2018) used VLT to measure the mass of Prox- 
ima Centauri through astrometric microlensing. Lu et 
al. (2016) employed the Keck telescope to look specifi- 
cally for isolated BHs by monitoring three microlensing 
events, where they achieved a final positional error of 
0.26 to 0.68 mas. The timescales of those events were 
60 to 160 days, and there were no detections of astro- 
metric deflections. 


3. IN SEARCH OF ISOLATED BLACK HOLES 
WITH AST 


3.1. Astrometry of Long-Duration Microlensing Events 


In 2009, we began a multi-cycle HST program of as- 
trometry of long-duration microlensing events in the di- 
rection of the Galactic bulge in order to detect isolated 
BHs and measure their masses. Our aim is to select 
events having timescales >200 days, light curves show- 
ing no evidence for a light contribution by a luminous 
lens, and preferably a high magnification factor. We 
then obtain high-resolution HST imaging as the events 
proceed, in order to measure the astrometric deflections 
of the background sources. To date we have monitored 
eight long-duration events. For some of them, there is 
no clear detection of an astrometric signal, but our data 
analysis is still in progress, and the results will be dis- 
cussed in separate publications. In the present paper 
we analyze and discuss our findings for an event that 
clearly shows a large astrometric deflection, consistent 
with a high-mass lens. 


3.2. MOA-2011-BLG-191/OGLE-2011-BLG-0462 


MOA-2011-BLG-191/OGLE-2011-BLG-0462  (here- 
after designated MOA-11-191/OGLE-11-462) was a 
long-duration and high-magnification microlensing 
event in the direction of the Galactic bulge. It was 
discovered independently by both MOA and OGLE 
ground-based microlensing survey programs, and an- 
nounced by both teams nearly simultaneously on 2011 
June 2, through their public-alert websites.! The tar- 
get was also covered by the Wise Microlensing Survey. 
Table 1 gives details of this remarkable event. 

MOA-11-191/OGLE-11-462 occurred in an extremely 
crowded Galactic bulge field, less than 2° from the 
Galactic center. The observed peak magnification factor 
of this event was only about 20 in the ground-based data, 
but this was strongly diluted by blending with neighbor- 


1 MOA alerts: — https://www.massey.ac.nz/~iabond/moa/ 
alerts. OGLE alerts: http://ogle.astrouw.edu.pl/ogle4/ews/ 
2011/ews.html 
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Table 1. Basic Data for MOA-11-191/OGLE-11-462 Microlensing Event 








Parameter 





Event designation (MOA) 

Event designation (OGLE) 

J2000 right ascension, a 

J2000 declination, 6 

Galactic coordinates, (J, b) 
Baseline F606W magnitude 
Baseline F814W magnitude 
Baseline (F606W — F814W) color 
Peak magnification, Amax 

Date of peak magnification, to 


Timescale, tz 








Value Sources & Notes* 
MOA-2011-BLG-191 (1 
OGLE-2011-BLG-0462 (1 
17:51:40.2082 (2 
—29:53:26.502 (2 
359° 86, —1°62 (2 
21.946 + 0.014 (3 
19.581 + 0.012 (3 
2.365 + 0.026 (3 
369 (4 
2011 July 20.825 (4 
270.7 + 11.2 days (4 











“Sources and notes: (1) MOA and OGLE websites; the event was first alerted by MOA; (2) This paper, from astrometric 
analysis in §5.2 in Gata EDR3 frame at average epoch 2013.5; (3) This paper, Vegamag scale, from photometric analysis in 


§5.3; (4) This paper, from Table 6. 


ing stars. It soon became apparent, based on findings 
disseminated through internal communications in the 
microlensing groups, that the undiluted event actually 
had an extremely high magnification factor, approaching 
400. Blending also made the apparent timescale of the 
event appear shorter than the actual value, which was 
inferred to be longer than 200 days. It was clear from 
the ground-based observations that there was blending, 
for two reasons. First, the light curve for a typical event 
has a characteristic shape that is completely determined 
by the timescale and the maximum magnification, ex- 
cept for distortions due, e.g., to the lens-source relative 
parallax. The shape of the observed light curve was in- 
consistent with the expected shape unless the light at 
baseline was highly diluted by a blend, thus implying 
that the real magnification was much larger than the 
observed value. Second, as the source brightened, its 
centroid position in the ground-based images was seen 
to change, again consistent with blending with a neigh- 
boring star. Note that this shift is due simply to blend- 
ing and scales with the separation of the two stars; it is 
unrelated to the much smaller relativistic deflection of 
the source itself, which is discussed below. 

Figure 4 shows an 8” x 8” region centered on the 
source, as imaged by us in the F814W (J-band) fil- 
ter by HST with its Wide Field Camera 3 (WFC3). 
The source star is encircled in green. A conspicuous 
neighbor, nearly 20 times brighter than the unmagni- 
fied source, lies at a separation of only 04. The cyan 
circles in the figure have diameters of 1” and 2”, cor- 
responding to the generally best seeing in the ground- 
based survey observations, and more typical seeing, re- 


spectively. Thus, in ground-based images, the source is 
indeed blended with the bright neighbor and a number 
of fainter stars, depending on the seeing. 

High-magnification events are generally very sensitive 
to perturbations due to planets around the lensing ob- 
jects (Mao & Paczyniski 1991; Griest & Safizadeh 1998). 
Thus considerable interest was aroused by MOA-11-191/ 
OGLE-11-462 among groups engaged in searches for 
such planets. As a result, intensive photometric mon- 
itoring of this event was carried out by multiple groups, 
providing valuable data for our analysis. 


4. HST OBSERVATIONS 


The MOA-11-191/OGLE-11-462 event satisfies all the 
selection criteria for our HST follow-up program de- 
scribed in §3.1, and thus we triggered our observing 
sequence. Our project had a “non-disruptive” target-of- 
opportunity status, requiring a lead time of about two 
to three weeks from activation to the first observations. 
The first-epoch HST data were obtained on 2011 Au- 
gust 8, some 19 days after the peak light magnification 
on 2011 July 20. The magnification was still reasonably 
high (~12, corresponding to u ~ 0.08), so that the ex- 
pected astrometric deflection was 6 ~ 0.04 0g [see Equa- 
tion (6)], i.e., close to zero at this epoch, but its correct 
value is taken into account in the model described in §8. 
Subsequent HST observations indicated departure from 
a linear proper motion for the source. Thus we con- 
tinued the imaging, ultimately over an interval of over 
six years, long enough for robust separation of the rela- 
tivistic deflection from proper motion. Table 2 gives the 
HST observing log. 
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Table 2. Journal of HST Wide Field Camera 3 Observations 











Epoch Date MJD Year Proposal No. Frames No. Frames 
ID in F606W* in F814W* 

1 2011 Aug 8 55781.7 2011.600 GO-12322 4 5 

2 2011 Oct 31 55865.2 2011.829 GO-12670 3 4 

3 2012 Sep 9 56179.2 2012.689 GO-12670 3 4 

4 2012 Sep 25 =56195.3  2012.733 + GO-12986 3 4 

5 2013 May 13 56425.8 2013.364 GO-12986 3 4 

6 2013 Oct 22 56587.2 2013.806 GO-13458 3 4 

7 2014 Oct 26 56956.1 2014.816 GO-13458 3 4 

8 2017 Aug 29 57994.7 2017.660 GO-14783 3 4 





“Individual exposure times ranged from a minimum of 60s at Epoch 1, to a maximum 


of 285 s at later epochs. 





Figure 4. HST image in the F814W (J-band) filter of an 
8” x 8” region centered on MOA-11-191/OGLE-11-462, ob- 
tained at our final epoch in 2017 August. North is at the top, 
east on the left. Encircled in green is the source star, now 
returned to baseline luminosity. The site is resolved into the 
source, a much brighter neighboring star 0°’4 to the WNW, 
and several nearby fainter stars. The inner cyan circle has 
a diameter of 1”, corresponding to the typical best seeing 
in ground-based microlensing survey images; the outer cyan 
circle’s diameter is 2”, which is not unusual seeing. The 
source, bright neighbor, and several fainter stars are gener- 
ally blended in ground-based frames, and the blending in- 
creases with seeing. 


All our HST observations were obtained with the 
UVIS channel of WFC3, whose CCD detectors provide 


a plate scale of 39.6 mas pixel—!. To avoid buffer dumps 
during the orbital visibility period and thus maximize 
observing efficiency, we used the UVIS2-2K2C-SUB sub- 
array, giving a field of view (FOV) of 80” x 80”. This 
FOV is large enough to provide dozens of nearby astro- 
metric reference stars surrounding the primary target. 

The WFC3 detectors are subject to an increasing 
amount of degradation of their charge-transfer efficiency 
(CTE) as they are exposed to the space environment. 
The chosen subarray aperture places the target in the 
middle of the left half of the UVIS2 CCD, which lessens 
the impact of imperfect CTE relative to a placement 
closer to the center of the FOV. Nevertheless, a time- 
dependent correction for CTE must still be applied in 
the astrometric analysis of the images. 

Our HST observations were taken at a total of eight 
epochs, strategically scheduled for measurement and 
characterization of the astrometric deflections. At each 
epoch, we obtained images in two filters (to verify the 
achromatic nature of the event, and to test for blend- 
ing by very close companions): “V” (F606W) and “J” 
(F814W). At the initial epoch, when the source was 
bright, we obtained nine exposures, four in F606W and 
five in F814W. At each subsequent epoch, using longer 
integration times because of the fading of the source, we 
obtained seven exposures, three in F606W and four in 
F814W. Individual exposure times were adjusted to take 
into account the brightness of the source and the orbital 
visibility of HST, and ranged from a minimum of 60 s 
at the first epoch to a maximum of 285 s at the later 
epochs. The telescope pointing was dithered by ~200 
pixels (~8”) between individual exposures; this allowed 
retention of a common set of reference stars in all the 
exposures, in order to mitigate errors in the distortion 
solution. To maximize the S/N for the most crucial as- 
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trometric measurements, we separated Epochs 3 and 4 
by only 16 days in 2012 September, around the time 
when the deflection was expected to be near maximum. 

Figure 5 zooms in on the field around the source in 
Figure 4, showing a 21 x 20 region as observed at all 
eight epochs. The bright source is marked with an ar- 
row in the Epoch 1 (top left) image, and it can be seen 
to fade in the subsequent frames. The astrometric de- 
flection was highest at Epochs 3 and 4, even though 
the photometric magnification was only about 10% at 
this epoch. There was very little photometric change in 
the subsequent epochs, but the astrometric deflections 
remained detectable until Epoch 7, demonstrating the 
need for astrometric monitoring over a much longer du- 
ration than the photometric-variability period. 


5. HST DATA ANALYSIS 
5.1. Image Processing 


We used the flat-fielded and CTE-corrected (_f1c) im- 
ages produced by the Space Telescope Science Institute 
pipeline reductions (Sahu et al. 2021; Dressel et al. 2021) 
for the analysis. As noted above, WFC3 suffered from 
increasingly poor CTE during this period, so it was 
essential to take it into account. The _flc products 
were produced using the v2.0 pixel-based CTE model 
described by Anderson (2021). 


5.2. Astrometric Analysis 


To measure stellar positions in individual frames, we 
used an updated version of the star-measuring algorithm 
described in Anderson & King (2006). The routine goes 
through each exposure pixel by pixel, and identifies as 
a potential star any local maximum that is sufficiently 
bright and isolated. The routine uses the spatially vari- 
able effective point-spread functions (PSFs) provided at 
the WFC3/UVIS website? to fit the PSF to the star 
images in the individual _f1c exposures, in order to de- 
termine a position and flux for each star in the raw pixel 
frame of that exposure. Finally, the positions are cor- 
rected for geometric distortion using the distortion so- 
lutions provided by Bellini et al. (2011). 

As the positions of individual stars are expected to 
change during the ~6-year course of our observations 
due to their proper motions, we needed to determine 
their proper motions to properly specify the reference 
frame. For this, we began with the Gaia Early Data 
Release 3 (EDR3) (Gaia Collaboration et al. 2021) po- 
sitions and motions for the bright but unsaturated HST 
stars in the field. The reference frame was constructed 


7 https://www.stsci.edu/hst/instrumentation/wfc3/ 


data-analysis/psf 


to place the bright star close to MOA-11-191/OGLE- 
11-462 at the center of the reference frame at (x,y) = 
(1000, 1000) at the 2016.0 epoch, with a plate scale of 
40 mas pixel”! and north up. (Note that the Gaia cat- 
alog could be incomplete in this region because of the 
high source density.) Using the Gaia positions and mo- 
tions, we determined the position for each Gaza star in 
this frame at each epoch in order to properly transform 
the distortion-corrected observations at that epoch into 
the reference frame. This ensures that the proper mo- 
tions that we derive represent absolute proper motions. 
After this initial set-up of the reference frame based on 
the brighter stars, we incorporated high-precision HST 
stars that were too faint to be found with high precision 
in the Gata catalog and solved for their accurate posi- 
tions and motions. We then used their time-dependent 
positions to improve the reference frame. Even after al- 
lowing individual solutions to improve on the basis of 
HST observations, there remains very good agreement 
between our proper motions and those of Gaia. Fig- 
ure 6 plots the proper motions of our reference stars 
derived from our HST observations against the Gaza 
proper motions, where the red points are for brighter 
stars with G < 18, and blue points are for fainter stars 
with G > 18. The agreement is imperfect, of course, 
since the HST observations have a 6-year baseline and 
have higher S/N in individual measurements, resulting 
in higher accuracy in proper-motion measurements, par- 
ticularly at fainter magnitudes. The agreement is better 
for the brighter sample, for which Gaia proper-motion 
errors are typically smaller (<0.2masyr~'). Note that 
the MOA-11-191/OGLE-11-462 source itself is too faint 
at baseline for inclusion in the Gaza catalog. 

In the next step of our analysis, we used only the 
HST observations, because the Gaia measurements have 
much higher uncertainties for the fainter stars, and the 
HST observations have a longer baseline of 6 years com- 
pared to the 3 years of Gaia. In this step of the trans- 
formation, we used stars (1) with brightness similar to 
the average brightness of the target, (2) with color sim- 
ilar to the source’s color (see Figure 7), and (3) lying 
within 350 pixels of the source. The first criterion min- 
imizes any residual shift caused by CTE effects. We 
note that we already used the most recent CTE correc- 
tion software for our analysis. Since the CTE effects 
on the position measurements are differential, using lin- 
ear transformations based on stars of similar brightness 
should remove any residual CTE effects. (It is worth 
noting here that the images with the highest astromet- 
ric deflection were taken when WFC3/UIVS was young, 
and when CTE losses were small.) The second criterion 
ensures that the stars used in the transformation belong 
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Figure 5. 2’'1 x 20 cutouts around MOA-11-191/OGLE-11-462 as observed by HST at all 8 epochs. Exposures were taken 
in F606W and F814W from 2011 to 2017; see Table 2 for details. The source star is marked by an arrow in the first-epoch 
(top left) image. At this epoch, on 2011 August 8, the magnification was a factor of ~12. The maximum astrometric deflection 
occurred at Epochs 3 and 4, when the photometric magnification was only about 10%. In 2017, the source had returned very 


close to its unmagnified brightness and undeflected position. 
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Figure 6. Proper motions of the reference stars used in 
our analysis derived from our HST observations with a 6- 
year baseline, versus the Gaia proper motions. Red points 
represent stars brighter than G = 18, for which the Gaia 
errors are typically <0.2 mas yr~', and blue points represent 
fainter stars for which the errors are larger. Our analysis 
uses the Gaia reference frame, so it is natural that there 
is good agreement between them. But the individual HST 
measurements have higher precision, particularly at fainter 
magnitudes. 


to the bulge, which helps in minimizing errors due to 
parallax, as described in more detail below. The third 
criterion minimizes residuals in the distortion solution. 

We employed an iterative procedure to measure the 
positions and proper motions of the stars, starting from 
the revised values in each iteration. We rejected the 
highest-sigma point after each complete iteration. We 


repeated this procedure until the highest-sigma point 
was no more than a preset tolerance, for which we 
adopted 60. Only a small number of points were re- 
jected by this procedure, mostly affected by cosmic-ray 
hits on the detector. Then at each epoch the reference- 
star positions were corrected for proper motion, and the 
positions of the source were determined relative to this 
adjusted frame. The estimated uncertainty in the posi- 
tion of the source star relative to the adjusted frame is 
~0.4 mas in each individual exposure. 

As an illustration, Figure 8 shows the proper mo- 
tions as measured for nine representative stars. Figure 9 
shows the errors in the proper motion measurements of 
the reference stars. We note that, all the reference stars 
are within about 1.5 mag of each other (See Figure 7). 
Figure 10 shows the the histograms of the residuals of 
each measurement from that star’s proper motion solu- 
tion along the RA and Dec directions. Both distribu- 
tions are consistent with a Gaussian distribution (the 
red curve). As shown by previous similar studies, the 
final reference-frame positions are expected to be inter- 
nally accurate to better than 0.01 pixel (Anderson et al. 
2008; Bellini et al. 2015). 

We specifically solved for the proper motions of the 
reference stars, but we ignored their parallaxes. The 
reason for adopting this approach is the following. 

Our choice of reference stars ensures that a large frac- 
tion of them belong to the Galactic bulge, and hence 
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Figure 7. Left panel: color-magnitude diagram [mreoew versus (Mreo6w — MFsi4aw)] for all stars in the HST/WFC3 field. The 
main-sequence turnoff occurs at MFe606w Z 22, above which the disk and the bulge split into two sequences: the redder stars 
are mainly bulge objects, while the bluer ones are mainly disk objects. Stars marked by blue squares are selected as “disk-like” 
stars, and the stars marked by black crosses are selected as “bulge-like.” The position of the (unmagnified) microlensed source 
is shown as a green dot. The nearby bright star 0'’4 away and about 20 times brighter than the unmagnified source is shown as a 
red dot. Magenta points are the astrometric reference stars used in our analysis. Top right panel: proper motions of the disk-like 
stars (blue), the reference stars (magenta), and the source (green dot). Bottom right panel: proper motions of the bulge-like 
stars (black crosses), the reference stars (magenta), and the source (green dot). Only bulge-like stars (shown in magenta) with 
brightness similar to the observed brightness of the star were used in the final astrometric transformations. This ensures that 


the reference stars are very similar to the target source star, and uncorrected CTE and small parallax effects should cancel in 
the differential astrometric measurements. 
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Motions of nine representative astrometric reference stars in right ascension and declination. The red lines are 


linear fits to the proper motions of the stars. 
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Figure 9. Proper-motion errors of the reference stars. 


The y-axis shows the proper-motion error, defined as o, = 
[ne + Tryl??, where oy2 and o,y are the proper-motion 
errors along the x and y axes, which are parallel to right 
ascension and declination, respectively. Reference stars cover 
a range of 1.5 mag around the (unamplified) magnitude of the 
source; their proper motion uncertainties vary significantly 
from object to object, with a modest systematic increase 
with magnitude. 
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Figure 10. Histogram of the position residuals for all the 
reference stars in RA and Dec from their proper motion solu- 
tion, as measured in each image separately, and scaled with 
the measured dispersion for that star. The scaled residuals 
closely follow a standard Normal distribution (shown by the 
red curves). As noted in the caption to Figure 9, reference 
stars cover a range of magnitudes around the unamplified 
source magnitude. 


have similar parallactic motion as the source star. We 
note that the source parallax is small to begin with (S0.2 
mas). Then the source position is referenced to stars 
chosen to be at comparable distance, so any remaining 
impact of the source parallax on the astrometry or pho- 
tometry of the event is expected to be negligible. 


As described earlier, there is a bright star ~10 pixels 
away from our source. The target star is close to the 
brightness of this neighbor in the first two epochs and 
slowly fades to its nominal brightness, at which point 
it is about 3 mag fainter than the bright neighbor. For 
accurate astrometry of the source, we wanted to make 
sure that the position measurements of the source are 
not affected by the presence of this bright neighbor. So 
we wanted to subtract the PSF of the bright star before 
measuring the positions of the source at every epoch. 
However, subtracting the bright star is not just a matter 
of subtracting a standard PSF. The separation is ~10 
pixels and the available library PSFs go out to 12 pixels, 
and are tapered and not very accurate in the wings. 
Thus we needed to make a more extended PSF model. 

To make an extended PSF, we carefully selected stars 
that (1) are within ~350 pixels of the bright neighbor, 
(2) have brightness and color similar to the neighbor, 
and (3) are fairly isolated. We found 18 such stars (ex- 
cluding the bright star itself), which provided a good 
sample to make the required extended PSFs. We used 
the images of these 18 stars to produce a separate well- 
sampled, extended PSF for each individual exposure. 
For illustration, Figure 11 shows the stacked PSFs in 
F606W and F814W. We have taken particular care to 
make sure that the PSF is well characterized in the wings 
since the source lies in the wings of the bright star, and 
subtracting the wings correctly is crucial for accurate 
astrometry. 

We then took this PSF model for each exposure and 
subtracted it from the neighbor star in each exposure. 
Figure 12 shows the original images (first and third 
rows) and the subtracted images (second and fourth 
rows) in the F606W and F814W filters. The residuals 
are very small, particularly in the wings of the PSF. We 
found that the astrometric position of the source changes 
by ~0.03 pixel (1.2 mas) after this subtraction, which 
could have a significant effect on the mass determination 
of the lens; so this extra step of neighbor subtraction was 
crucial in improving the analysis/results. The resulting 
astrometric positions of the source were used for further 
analysis as described in the next section. 


5.3. Photometric Analysis 


In addition to the measured positions, the analysis al- 
gorithm provides PSF-based photometry of all the stars 
in the field. To set a calibrated zero-point, we used 
standard aperture photometry to determine the fluxes 
of a few isolated stars in the field within an aperture 
with a 10-pixel radius. These fluxes were then corrected 
to an infinite aperture, using encircled-energy measure- 
ments from Calamida et al. (2021), and the photometric 
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Figure 11. The stacked F606W (left) and F814W (right) 
PSFs, with the 10-real-pixel radius shown in green. Note 
that the F606W PSF has a diffraction spike bump very close 
to the location of the source star. The bump for the other 
diffraction spikes does not show up in the residual images 
shown in the second row of Figure 12, which implies that 
it is subtracted well beneath the source, making the posi- 
tion measurements of the source more robust. The F814W 
PSF has a very strong radial gradient (along with azimuthal 
structure) at the location of the source. Without subtracting 
a high-fidelity model PSF, there would be some impact on 
the measurement of the source positions. And as the source 
moves relative to the neighbor, then the source would move 
across the PSF halo, which could introduce artificial shifts if 
the PSF of the bright neighbor is not subtracted well. 


zero-point in the image headers (PHOTFLAM) was used to 
convert these fluxes to the Vegamag scale. The mean 
difference between these values and the values obtained 
by the PSF fitting was then applied to all of the PSF 
magnitudes to convert them to Vegamag. 

As described above, the source star lies on the wings of 
the PSF of the neighboring bright star. So, for accurate 
photometry of the source, it was critical to correctly 
subtract the light contribution from the bright neigh- 
bor. The photometry for the source was carried out 
after subtracting the superposed flux from the neighbor 
star using a high-fidelity PSF as described above. The 
resultant time-series HST photometry of the source is 
shown in Figure 13, along with the model light curve de- 
scribed below in §8. There is no detectable color change 
as the event progresses and the star fades: the color of 
the source has remained constant to within 0.01 mag 
during the 6 years of observations with HST. There 
is also no detectable blending as described in more de- 
tail in 88.2. The source at baseline brightness has ap- 
parent magnitudes of mpeogw = 21.946 + 0.014 and 
Mrsiaw = 19.581 + 0.012. 

We made a stack of all the images for each filter for 
every epoch. The HST images allow us to detect and 
measure magnitudes of stars as faint as V ~ 25. Figure 7 











shows the CMD based on this photometry, where we also 
show the position of the source and the bright neighbor 
04 away. 


6. GROUND-BASED LIGHT CURVE 


MOA-11-191/OGLE-11-462 was monitored photomet- 
rically by several ground-based observatories. The cov- 
erage by MOA, OGLE, and Wise Microlensing Sur- 
vey extended over several years. Moreover, as a high- 
magnification event, it attracted intensive monitoring 
by a number of additional ground-based telescopes— 
especially around the time of peak brightness, where 
the microlensing light curve is sensitive to planet detec- 
tion. Table 3 gives a journal of the photometric obser- 
vations and data used in our analysis. We use the data 
re-reduced by the surveys and other groups (Udalski et 
al. 1992; Bond et al. 2001; Sackett et al. 2004; Gould et 
al. 2006; Tsapras et al. 2009; Dominik et al. 2010). 

Figure 14 shows the light curve, both over a 300-day 
interval (top panel), and zooming in on the seven days 
around peak magnification (bottom panel). Superposed 
is our model fit to the light curve, from the analysis 
described in the next two sections. 


7. BLENDING, RELATIVE PARALLAX, AND 
LENS TRAJECTORY 


The combination of ground-based photometric moni- 
toring with long-term astrometric and photometric mea- 
surements from HST affords the ability to constrain all 
aspects of this event and obtain high-quality measure- 
ments of its parameters. In this section, we describe 
some characteristics of the data, illustrating the key in- 
formation that can be obtained from photometry and 
astrometry separately through heuristic considerations. 


7.1. Photometric Blending 


Ground-based photometry of MOA-11-191/OGLE-11- 
462 suffers from significant amounts of blending with 
neighboring stars, as shown in our HST images (Fig- 
ure 4). Moreover, the amount of blending changes 
markedly depending on image quality. There is a bright 
neighbor star only 0’4 away, along with two fainter stars 
within 05 of the source. At larger separations, there are 
three more stars within 1”, whose combined brightness 
is greater than that of the baseline source, and there are 
several more stars within 15 which are also brighter 
than the source. Since in the available ground-based 
imaging data the measured image quality is seldom bet- 
ter than 1”, the ground-based photometry will always 
include the light from at least the three closest stars 
within 0/5. In a fraction of the data, taken under poorer 
seeing conditions, the source is blended with an increas- 
ing number of neighbors. 
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Figure 12. The source star lies in the PSF wings of a bright neighbor, marked in these frames with green circles with a radius 
of 10 pixels. Since the available WFC3 library PSFs do not extend to this large a radius, a special PSF extending to 20 pixels 
was constructed using the same HST images. This PSF was used to subtract the neighbor before measuring the position of the 
source in each exposure. The top row shows the original stacked F606W images of Epochs 1 to 8 (from left to right), and the 
second row shows the PSF-subtracted images. The third row shows the original stacked F814W images, and the bottom row 


shows the subtracted images. The stacks are overbinned by a factor of 2 to show details. Care was taken to assure that the 
PSF is well characterized in the wings of the neighbor star where the source lies. 


Our HST images allow us to place constraints on the 
expected blending parameter, g, defined as the ratio 
of the flux from neighbors included in the photometry 
to the flux from the unmagnified source itself [Equa- 
tion (5)], in the ground-based observations. The bright 
neighbor is 18.88 times brighter than the source at base- 
line in F814W, and the contribution from the two fainter 
stars is an additional 0.19 times that of the source. Thus 
the expected blending factor due to these three stars is 
g = 19.07 in the F814W filter. Since the bright neigh- 
bor is similar in color to the source (see Figure 7), we 
adopt this value of g for both OGLE (J band) and MOA 
(R band) data as an initial estimate, but keep it as a 
variable in our analysis. The final values (§5.2) differ 
significantly between OGLE and MOA, possibly due in 
part to differences in processing between the two data 
sets. 

The effect of blending can be reduced substantially 
by basing the photometry on difference images. It can 
be further reduced by restricting the analysis to images 
taken under good seeing. It is obvious, however, that 
variable blending will affect the noise characteristics of 
ground-based photometry; it is difficult to include such 
effects deterministically because of the imperfect knowl- 


egde of the blending (at the sub-percent level) for indi- 
vidual images. 


7.2. Heuristic Considerations 


7.2.1. Photometric Constraints on Parallax and Lens 
Trajectory 


As discussed in §2.3, the light curve of a long-duration 
microlensing event such as MOA-11-191/OGLE-11-462 
can show distortion by the relative parallactic motions 
of the source and lens (e.g., Gould 1992; Alcock et al. 
1995). Specifically, the light curve is sensitive to 7 
and yrs, because their combination modifies the relative 
path of source and lens, and thus the shape of the light 
curve. We note that, in our formalism, yig corresponds 
to the position angle (PA) of the path of the lens relative 
to the source without parallax in equatorial coordinates 
(not to be confused with the instantaneous path of the 
lens at the time of closest angular approach in ecliptic 
coordinates; see §2.2). 

In principle, a sufficiently accurate light curve can pro- 
vide good constraints on both mg and yrs. However, 
as discussed in the previous subsection, the photometry 
is significantly affected by blending. We attempted to 
model the light curve alone, but found that it can be 
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Table 3. Journal of Ground-Based Photometry of MOA-11-191/OGLE-11-462 











Data Set Telescope Aperture Filter No. of Date Range 
Location [m] Observations (HJD — 2450000) 
MOA New Zealand 1.8 R 46040 3824.093...9441.117 
OGLE Chile 1.3 I 15546 5260.855. . .8787.509 
Wise Survey Israel 1.0 I 953 5658.529...5722.543 
Danish DFOSC Chile 1.54 I 921 5744.804...5782.554 
Danish LuckyCam Chile 1.54 broad 10 5738.624...5742.706 
MONET North Texas, USA 1.2 I 214 5762.722...5764.824 
Faulkes North Hawaii, USA 2.0 SDSS 3’ 99 5763.777 ...5768.955 
Liverpool Canary Islands, Spain 2.0 SDSS 3’ 254 5739.521...5768.438 
SAAO 1.0m South Africa 1.0 I 611 5751.264...5777.315 
SAAO 1.0m South Africa 1.0 V 44 5758.454...5765.385 
U. Tasmania Australia 1.0 I 60 5761.055...5769.110 
CTIO Chile 1.3 I 226 5757.511...5772.785 
CTIO Chile 1.3 V 27 5757.515.. .5763.807 
Auckland New Zealand 0.4 R 160 5759.873 ...5778.893 
Farm Cove New Zealand 0.35 unfiltered 37 5741.856...5761.856 
Kumeu Obs. New Zealand 0.35 R 63 5759.803...5762.134 
Vintage Lane New Zealand 0.4 unfiltered 60 5762.792...5767.941 
Weizmann Israel 0.4 I 167 5762.281...5764.423 
Wise Israel 0.46 I 142 5762.266...5763.456 





Note—The data were acquired by the MOA (Bond et al. 2001), OGLE (Udalski et al. 2015), Wise Microlensing Survey 
(Shvartzvald et al. 2016), MiNDSTEp (Dominik et al. 2010), RoboNet (Tsapras et al. 2009), PLANET (Sackett et al. 2004), 


and «FUN (Gould et al. 2006) teams. 


fitted with a range of parameter combinations, in which 
the values of mg and yzs are strongly correlated. In ad- 
dition, the derived value of 7_ varies with the specific 
subset of photometric data chosen for analysis, as well 
as with the assumed blending factor for those data. The 
derived value of 7g ranges from 0.07 to 0.12, with larger 
values corresponding to larger values of yrs, ranging 
from 330° to 358°. The reason is that increasing ys 
makes the lens move in a more northerly direction as 
seen in the bottom panel of Figure 15. Since parallax 
is predominantly in the east-west direction, in order to 
produce a fixed change in u, the value of 7g has to in- 
crease with ys, so that the change in position due to 
parallax can compensate for a more northerly motion of 
the lens. Several different combinations of these quanti- 
ties can reproduce the observed light curve, with differ- 
ences between solutions of the order of 1 mmag at early 
and late times, and ~5 mmag near the peak. System- 
atic differences in the data at this level could be caused 
by small variations in blending associated with changes 
in the ground-based seeing, or other minor secular vari- 
ations in the photometry. Therefore we conclude that, 


when photometry alone is used to constrain the param- 
eters of the event, only a reliable joint constraint on 7p 
and yrs can be derived. Fortunately, astrometry pro- 
vides a robust independent estimate of yrs, allowing us 
to break this degeneracy and determine the two quanti- 
ties separately. 


7.2.2. Astrometric Deflection and Orientation of the 
Relative Motion 


In order to understand how astrometry can constrain 
the direction of motion of the lens, it is useful to consider 
an illustrative plot of the motion of the lens relative to 
the source, as shown in Figure 15. North is at the top, 
east on the left; the (uRA, UDec) coordinates give the po- 
sition of the lens relative to the source in units of 0p, 
with wra increasing to the east. We have used our ac- 
tual final model described below in §8.2 for this illustra- 
tion. The top panel shows the motion of the lens, with 
Te = 0.0894 and yrs = 342°5, and an impact param- 
eter of uo = 0.00271. The straight line represents the 
proper motion of the lens with respect to the source, 
while the wavy line adds the parallactic motion, com- 
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Figure 13. Top panel: Photometry of MOA-11-191/OGLE- 
11-462, obtained with HST over 6 years at 8 epochs in the 
F606W (V) and F814W (J) filters (black filled circles), along 
with our final model fits from §8.2 (red and blue curves). Er- 
rors on the photometry are smaller than the plotting points. 
The second panel plots the observed values of V —J, showing 
that the source color remained constant to within ~0.01 mag 
during the entire 6-year duration. The third and fourth pan- 
els show that the residuals, Vops — Vmoder and Ions — Imoael, 
at the epochs of the HST observations are consistent with 
zero within the measurement uncertainties (see §8.2 for more 
details and discussion). 


puted using the JPL ephemeris of the Earth.? Red dots 
show the position of the lens at the eight epochs of our 
HST observations. 

The bottom panel in Figure 15 shows an enlarged view 
of the lens trajectory near the source position. The 
dotted black line represents the proper motion of the 
lens with respect to the source without parallax, while 
the solid black line includes the parallactic effect. (The 
red lines correspond to a less-preferred uo,, solution de- 
scribed in §8.3.) The plot shows that near the closest 
angular approach—and thus the peak magnification— 
the relative path is substantially affected by parallax; 
however, the astrometric deflection is very small at this 
time (see §2.3). Since the source deflection is always in 
the direction of the line joining the instantaneous posi- 
tion of the lens to the undeflected position of the source, 
the directions of the source deflections at late times will 
remain nearly constant with little parallax effect; thus 


3 https: //ssd.jpl.nasa.gov /horizons/app. html / 


the late-time deflection directions robustly constrain the 
orientation, yrs, of the lens trajectory. 


7.2.3. Constraining the Lens Trajectory Orientation 


Since the parallactic effect is unimportant for con- 
straining yrs, we first fitted the photometry using a 
light-curve model that neglects parallax. The resultant 
model (with to — 2450000 = 5763.33, tm = 231.56 days, 
Amax = 372.62, and g = 19.3) predicts the total de- 
flection in units of Op as a function of time, through 
Equation (6). 

As described in §5.2, we have accurate measurements 
of the (x,y) positions of the source at the eight epochs 
of HST observations. These positions are affected both 
by the proper motion of the source, and its deflections, 
which are a function of Og and yrs. We fitted a model 
to the positions, whose parameters are the x (RA) and 
y (Dec) components of the proper motion, @n, and yrs. 
This fit resulted in values of Og = 5.2 + 0.5 mas and 
pis = 337°9 + 5°0. 

Components of the resultant deflections as a function 
of time, after subtracting the best-fitting proper motion, 
are shown in the top two panels of Figure 16. These pan- 
els plot the deflections in the RA and Dec directions. 
The bottom-left panel shows the total amount of deflec- 
tion, again as a function of time. Note that the total 
deflection reaches a maximum of ~2 mas in late 2012. 

The bottom-right panel of Figure 16 plots the RA 
versus Dec deflections. These deflections are always 
along the line joining the lens to the source, and thus at 
large separations their direction is opposite to the direc- 
tion of the relative motion of the lens. The solid black 
line passes through the origin, at the pig angle derived 
above. The dotted blue lines indicate the allowed range 
based on the uncertainties. Although the best fit value 
from photometry alone is yyy = 354°8 (shown by the 
solid magenta line), the allowed range of yrs (shown by 
the dashed magenta lines) has a flat probability distri- 
bution. The combined constraint from photometry and 
astrometry is used in our subsequent analysis. 

Note that for clarity in Figure 16 we have shown the 
mean deflections at each epoch, and not the individual 
measurements. In particular, this avoids a confusing 
overlap of points in the bottom-right panel, where the 
deflections are not a monotonic function of time. The 
individual measurements are shown in the next section. 


8. FULL MODELING OF THE PHOTOMETRIC 
AND ASTROMETRIC DATA 


In this section we give full details of our analysis, car- 
ried out independently by several coauthors using dif- 
ferent parameterizations, all leading to a consistent final 
model of the event. 
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Figure 14. Ground-based photometric observations of MOA-11-191/OGLE-11-462, along with a best-fitting model light curve, 
shown over a 300-day interval in the upper panel and over a zoomed-in region covering seven days around peak magnification in 
the bottom panel. All data have been transformed to OGLE J magnitudes according to the inferred baseline magnitudes and 
blend ratio from the common model. 


8.1. First Approach: All Photometric Data Sets, and icant constraints on the lens-source model—especially 
Robustness of Parallax Measurement since each set of observations can have a different base- 
In addition to OGLE and MOA, photometric time- line magnitude and blending parameter. Nevertheless, 


in our initial approach, we included all the data in our 
analysis, but also carried out analyses separately for 
“OQGLE-only,” “MOA-only,” and “OGLE+MOA-only” 
data sets. 

Ground-based time-series photometry is susceptible to 
systematic noise, and this must not be mistaken for real 


series data were obtained at several more observatories 
(see Table 2). These data typically cover a relatively 
narrow range of ~20 days around the peak, being pri- 
marily aimed at searching for planet-related distortions 
of the light curve. These photometric series, unlike those 
obtained by the survey programs, do not provide signif- 
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Figure 15. Top panel: path of the lens with respect to the 
source, the position of the source being fixed at (0,0). The 
direction of motion is shown by an arrow. The black straight 
line is the lens path without parallax, and the wavy line is 
the path including parallax. Red points mark the epochs of 
HST observations. urna and upec are the RA and Dec com- 
ponents of the lens-source separation, in units of 0g. Bottom 
panel: enlarged view of the lens trajectory with respect to 
the source around closest angular approach. Black lines cor- 
respond to the uo,— solution, which is the preferred solution, 
and the red lines show the less-preferred uo,+ solution. In 
both cases, the barycentric trajectories (no parallax), shown 
by dotted lines, pass on the north side of the source. The 
geocentric trajectory (i.e., the trajectory as seen by an ob- 
server on Earth, which includes parallax) of the uo,+ solution 
(shown by the solid red line) passes on the north side of the 
source, but the trajectory of the uo,— solution (shown by the 
black solid line) passes to the south. The wo,— solution is 
used in our analysis presented here. It is worth noting, how- 
ever, that the derived 6g, and hence the mass of the lens, is 
nearly identical for both solutions. yrs is the position angle 
of the barycentric proper motion of the lens with respect to 
the source (i.e., no parallax), as shown here. 


features of the light curve. To improve the robustness 
of our solution, we model the photometric uncertainties, 
and moreover force the model to follow the bulk of the 
data by explicitly down-weighting outliers. As imple- 
mented for the SIGNALMEN microlensing anomaly de- 
tector (Dominik et al. 2007, 2019), we specifically adopt 
a bi-square weight function with regard to the median 
residual, and a Gaussian distribution for the uncertain- 
ties with revised standard deviation in magnitude of 


6 = (no)? +03, (17) 


where o denotes the reported error bar, « is a scaling 
factor, and op corresponds to a systematic error added 
in quadrature. For the plot of the various data sets as 
shown in Figure 14, we give the respective estimated 
values of & and og in Table 5. If the size of the error 
bars does not vary substantially, there is a degeneracy 
between « and go, and either of the parameters provides 
modified constant error bars (while it does not matter 
which). 

We find two viable models, significantly only distin- 
guished by the sign of ug, and will refer in the following 
to the model with uo,_. The microlensing parallax pa- 
rameter 7p is constrained by the wing of the light curve 
and much less sensitive to the peak region, and therefore 
is mostly constrained by the microlensing survey data. 
From various combinations of data sets, we consistently 
find mz = 0.10 + 0.02, but some variation in the trajec- 
tory angle w, correlated with tg and the blend fraction, 
yielding visually indistinguishable model light curves. 

However, the angle of lens-source proper motion, yrs, 
follows robustly from the astrometric data (see §7.2.2), 
given that the centroid shift to first order (i.e., neglect- 
ing the small distortion caused by parallax) traces an 
ellipse (a highly flattened ellipse resembling a line in 
our case) whose semi-major axis is parallel to ps. If 
we restrict this angle to the range 333° < yrs < 348°, 
as suggested by the astrometric data, the photometric 
light curve does not change substantially, and we find 
TE & 0.086. 

We emphasize here that it is incorrect to say that 
there is a discrepancy between the paths determined 
from photometry and astrometry, since there is a cor- 
relation between yrs and 7p in the photometric solu- 
tion. However, restricting the trajectory angle ys as 
robustly derived in the last section from the orientation 
of the centroid shifts to the range 333° < yrs < 348°, 
suggests 7p 0.086. 





8.2. Second Approach: Simultaneous Fit of 
Photometric and Astrometric Data 
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Figure 16. The top two panels show the average values of the measured deflections in RA and Dec at each HST epoch. We 
use the values of to,tm, and uo as derived from a light-curve fit without parallax, and fit for the proper motion of the source, 
Oz, and yys. The solid black line is the best fit with 0g = 5.2 mas and yr_s = 337°9. The bottom left panel shows the total 
deflections at each epoch. The bottom right panel shows the RA versus Dec deflections. These deflections are always along 
the line joining the lens to the source, and thus at large separations their direction is opposite to the direction of the relative 
motion of the lens with respect to the source. The solid black line passes through the origin, at the yis angle derived above. 
The dotted blue lines indicate the allowed range from astrometry based on the uncertainties. The allowed range of yis from 
photometry alone is shown by the dotted magenta lines, the best fit value being yrs = 354°8 shown by the solid magenta line. 
The combined constraint from photometry and astrometry is used in our subsequent analysis. 
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We now turn to a full analysis in which we fit the as- 
trometric and photometric data simultaneously in order 
to obtain all of the parameters. Such a solution is im- 
portant, since the crucial parameters of 6g and 7p are 
derived from two different types of data. A simultane- 
ous solution is also essential for a correct estimate of the 
uncertainties in the model parameters. 

We follow the same plane-of-the-sky approach de- 
scribed in §7.2, which makes it easier to work with, and 
also show the actual paths of the lens and the source, 
and the deflections. We follow a different parameteriza- 
tion procedure where the model parameters we optimize 
contain all terms needed to characterize the positions 
of the lens and the source on the sky as a function of 
time; these include the reference positions and proper 
motions of both lens and source, their relative parallax, 
and the angular Einstein radius of the lens. In prin- 
ciple, the source parallax is also needed; however, its 
parallax in the reference system we use is close to zero, 
and it is not meaningfully constrained by the observa- 
tions. As discussed below in §9.2, the best constraints 
on the source distance come instead from photometry 
and high-resolution spectroscopy. 

From these parameters, the undeflected paths of the 
source and lens can be determined. The deflection of 
the lensed image of the source is then computed, and 
the resulting deflected source positions are matched to 
the observed positions. The same calculation also yields 
the source magnification; in order to match the observed 
photometry, the model must include a baseline magni- 
tude of the source and a blending parameter for each 
photometric data set. Consistent with the previous ap- 
proach, we found that most of the photometric data 
sets cover too short a time interval to yield meaning- 
ful constraints on the event parameters in the presence 
of significant blending; therefore we limit the model op- 
timization to the MOA and OGLE photometric data 
sets, and validate the resulting model for the other data 
sets separately (see §8.1). Also, in order to avoid un- 
due impact from any secular variations in photometric 
responses, we only include OGLE and MOA photomet- 
ric measurements within +2 years from the peak of the 
event. We adopt the approximate values of yzs, tp, 
to, and ug from the analysis of the previous section as 
our initial estimates, but we leave all parameters free 
in the optimization; the baseline magnitude and blend- 
ing parameters for MOA and OGLE are also separately 
optimized. 

As discussed in the previous subsection, the results of 
the optimization depend to some extent on the relative 
weighting of astrometry and photometry. Because the 
number of photometric measurements greatly exceeds 





that of astrometric measurements, and the nominal pho- 
tometric uncertainties are very small, an optimization 
using nominal errors disproportionately weights pho- 
tometry, resulting in a poor match to the astrometry. 
In order to obtain a more balanced weighting of astrom- 
etry and photometry, we scaled the photometric errors 
by different amounts in different temporal bins, and val- 
idated each solution based on a reasonable match to the 
astrometric data. In our final model, most photomet- 
ric points are still at an uncertainty below 10 mmag. 
This solution has a total astrometric x? of 136 with 
106 points. (The solution with nominal weights has a 
higher astrometric 7 of 149.) The parameters of the 
final model are given in Table 4. 

Figure 17 shows the reconstructed motion of the lens 
(magenta) and of the source (black) in the plane of the 
sky, based on our adopted model. The predicted appar- 
ent source trajectory is shown by the green solid line. 
The astrometric measurements are shown individually 
by the small filled circles, and their epoch averages as 
red triangles. Cyan squares show the model position at 
each epoch; gray lines connect the model lens position 
to the undeflected and deflected source positions at each 
epoch. 

Figure 18 presents the measured and predicted source 
positions separately for the x and y coordinates. To im- 
prove the legibility of the plot, the fitted proper motion 
of the source has been subtracted from both model and 
data; therefore the points shown represent the deflection 
of the source. The black line is our final adopted model 
which takes photometry as well as astrometry into ac- 
count. 

We can now check the consistency of the final model 
with our measured HST photometry in the F606W and 
F814W filters, and also constrain the amount of blend- 
ing. We used the final model to calculate the magnifi- 
cations at the HST observation epochs, and varied the 
baseline magnitudes and the blending factor to fit the 
HST photometry. The resulting fit is excellent (see 
Figure 13), and yields baseline source magnitudes of 
Mpeoow = 21.946+40.012 and mpgiaw = 19.581 40.012, 
with corresponding blending factors for the HST pho- 
tometry of g = —0.012 + 0.015 and —0.006 + 0.012, re- 
spectively. The bottom two panels of Figure 13 show 
the residuals of the observations relative to the model 
in F606W and F814W at the epochs of HST observa- 
tions, where we have assumed the minimum physically 
allowed value of g = 0. (It makes little difference if we 
instead use the slightly negative values of g from the 
formal fit.) The stringent constraints on blending and 
lack of color variation make it unlikely that our deflec- 
tion measurements could be affected by blending with a 
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Table 4. Parameters of the Full Fit to Astrometry and Photometry 














Parameter Units Value Uncertainty (lo) Notes* 

pis (RA) mas yr + —2.263 0.029 (1) 
tus (Dec) mas yr + —3.597 0.030 1) 
On mas 5.18 0.51 (2) 
th days 270.7 11.2 3) 
LS deg 342.5 4.9 (4) 
t6 (HJD—2450000.0) days 5765.00 0.87 5) 
TE 0.0894 0.0135 (6) 
Uo 0.0422 0.0072 7) 
MOA Baseline R magnitude mag 16.5147 0.0016 

MOA Blending parameter 16.07 0.66 

OGLE Baseline J magnitude mag 16.4063 0.0015 

OGLE Blending parameter 18.80 0.79 

Derived parameters: 

uo 0.00271 (8) 
to (HJD—2450000.0) days 5763.32 (9) 





“Notes: (1) Undeflected proper motion of the source in Gaia EDR3 absolute frame; 
(2) Angular Einstein radius; (3) Angular Einstein radius 0g divided by absolute value 
of lens-source proper motion pits; (4) Orientation angle of the lens proper motion 
relative to the source (N through E); (5) Time of closest angular approach without 
parallax motion; (6) Relative parallax of lens and source in units of 6m; (7) Impact 
parameter in units of Og without parallax motion; (8) Impact parameter derived using 
the model-fit parameters above and after including parallax motion, in units of 05; 
(9) Time of closest angular approach derived as above and after including parallax 


motion. 


binary companion or a field star lying within the HST 
PSF. 


8.3. Lens Motion 


The path of the lens in the sky plane, as derived from 
the above analysis, is shown in Figure 15. We note that 
there are two solutions corresponding to uo,+ and wo,— 
(Gould 2004). However, at to, the angular separation 
between the source and the lens is dominated by the 
parallactic motion of the lens (the separation caused by 
the parallactic motion is ~0.04, compared to the actual 
value of ug Y 0.00271). Thus, the path without parallax 
in both solutions lies on the uo,+ side. 

The paths of the lens for the uo,+ and uo,— solutions 
are shown in Figure 15, where the dashed lines show 
the path of the lens without parallax, and the solid lines 
show the path of the lens after the parallactic motion is 
taken into account. The respective paths are separated 
only by 0.02 mas at the time of maximum magnification, 
and quickly merge. The uo,— solution is the preferred 
solution, and is the one used in our analysis. However, 
we verified that both solutions provide practically iden- 


tical deflections (since the deflection measurements are 
at much higher u, where the two paths nearly merge), 
and the results are the same for all practical purposes 
in both solutions. 

The first two rows of Table 4 give the proper motions 
of the source. We use the values of 05, t§, and yrs given 
in the next three rows to determine the proper motion 
of the lens with respect to the source as —2.10 + 0.22 
and 6.66 + 0.67 mas yr~! in RA and Dec, respectively. 
The resulting absolute proper motions of the lens are 
given below in Table 6. 








9. PROPERTIES OF THE SOURCE 


As described in §§2.2—2.3, the mass determination for 
the lens does not depend upon the individual distances 
to the lens and source, but only on the relative lens- 
source parallax, mpg, and the Einstein ring radius, 9g— 
quantities that are directly determined from the light 
curve and the measured astrometric deflections. How- 
ever, we still need an estimate of the distance to the 
source, Dg, in order to determine the distance to the 
lens, D,, which is discussed below. Moreover, MOA-11- 
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Figure 17. Representation of the reconstructed motions of the lens and of the source on the plane of the sky. The reference 
point is the undeflected position of the source at time to. The small dots are the individual measurements from HST images 
(green for F606W, orange for F814W)); the larger red triangles represent the average positions for each epoch. The cyan squares 
are the fitted positions at each epoch. The gray lines connect the undeflected source and lens positions at each HST epoch. 


191/OGLE-11-462 was a very high-magnification event, 
with an impact parameter of only uo ~ 0.00271. Thus 
it is desirable to estimate the angular diameter of the 
source, to verify that it is consistent with the point- 
source light-curve modeling adopted in the previous sec- 
tion. 

In this section we use results of ground-based spec- 
troscopy of the magnified source, and HST photometry 
at its baseline brightness, to estimate its distance and 
angular diameter. 


9.1. High-Resolution Spectroscopy 


In a target-of-opportunity program focusing on high- 
dispersion spectroscopy of high-amplitude microlensing 
events in the Galactic bulge, Bensby et al. (2013) ob- 
tained four spectra of MOA-11-191/OGLE-11-462 on 
2011 July 20 and 21. These observations were made 
almost precisely around the dates of maximum magni- 


fication. The authors employed spectrographs on three 
different large telescopes: VLT, Magellan, and Keck I, at 
spectral resolutions ranging from 46,000 to 90,000. The 
aim of their program was to determine chemical com- 
positions of dwarfs and subgiant stars in the Galactic 
bulge. 

Bensby et al. (2013) carried out a non-LTE model- 
atmosphere analysis of the spectra, obtaining an effec- 
tive temperature and surface gravity of Teg = 5382 4 
92 K and logg = 3.80 + 0.13, consistent with a late G- 
type subgiant. The metallicity was found to be slightly 
above solar, at [Fe/H] = +0.26 + 0.14, and the radial 
velocity is +134.0 kms7!. 

As discussed above (e.g., Figures 4 and 5), the HST 
images show that MOA-11-191/OGLE-11-462 is accom- 
panied by a neighboring star only 0/4 away, which would 
have been included in the spectrograph apertures. How- 
ever, at the dates of the observations, MOA-11-191/ 
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Figure 18. Predicted and measured positions of the source for our adopted joint astrometric and photometric fit. Small dots 
show the individual measurements from HS'T images (green for F606W, orange for F814W), while red triangles with error bars 
show the average and uncertainty at each HST epoch. The black line is our final adopted model. The fitted proper motion for 
the source has been subtracted from both model and measurements, in order to allow a better scaling of the plot. 


OGLE-11-462 was ~3.2 mag brighter than the neigh- 
bor in the F606W bandpass; thus the contamination 
was about 5%, which should have a minimal effect on 
the spectroscopic investigation. 

As shown in Figure 7, the neighbor has a color similar 
to that of the baseline source, but is more luminous. We 
verified its small effect on the spectroscopy by artificially 
contaminating the source spectrum with a 5% flux from 
a giant star with the same effective temperature and 
metallicity, but with log g = 0. This produced a smaller 
change in the line strengths of the combined light than 
would a change of log g for the source star by its 0.13 dex 
uncertainty. 


9.2. Distance and Angular Diameter 


We estimate the distance to the source using two inde- 
pendent sets of theoretical stellar isochrones: “Padova” 
isochrones, which we obtained using the PARSEC web 
tool* (Bressan et al. 2012; Marigo et al. 2017); and 
“BaSTI” isochrones? (Hidalgo et al. 2018; Pietrinferni et 
al. 2021). In Figure 19 we plot the position of the source 
in the distance-independent logg versus log Teg plane 
(red filled circle with error bars). Superposed are PAR- 
SEC isochrones, calculated for the measured metallicity 


4 Version 3.6 at http://stev.oapd.inaf.it/cmd 
5 http: //basti-iac.oa-abruzzo.inaf.it /index.html 
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Table 5. 
Data 


Revised Error Bars of Photometric 








Data set 





OGLE 
MOA 
Wise Survey 


Danish 1.54m DFOSC 
Danish 1.54m LuckyCam 


MONET North 1.2m 
Faulkes North 2.0m 
Liverpool 2.0m 
SAAO 1.0m I 
SAAO 1.0m V 
Tasmania 1.0m 
CTIO 1.3m I 
CTIO 1.3m V 
Auckland 0.4m 
Farm Cove 0.35m 
Kumeu 0.35m 
Vintage Lane 0.4m 
Weizmann 0.4m 
Wise 0.46m 


K 00 
1.25 0.005 
1.18 0.0025 
0.64 0.009 
2.75  10°-° (*) 

0.1 (*) 0.050 
0.57 0.006 
0.1 (*) 0.014 
1.17 0.008 
1.42 0.0014 
0.59 0.009 
0.28 0.017 


42 107° @) 
1.01 107° (*) 
0.91 0.007 
1.00 0.003 
1.87 107° (*) 
Cie; wos 
0.80 107° (*) 
0.22 0.009 





NotE—The adopted error bar (in magnitude) becomes & = \/(K a)? + 02, where o is the reported error bar of the photometric 
data. We have applied range constraints « > 0.1 and oo > 107°. An asterisk (*) indicates that the value is at the range 


boundary. 


of [M/H] = +0.26 (blue curves), and BaSTI isochrones, 
for a metallicity of [M/H] = +0.3 (orange curves). Both 
sets of isochrones agree well around the location of the 
source star, yielding an age of about 4 Gyr, but with 
considerable uncertainty. 

At the position of the source in the diagram, a star on 
the PARSEC 4-Gyr isochrone has absolute magnitudes 
(Vegamag scale) in the native HST/WFC3 bandpasses 
of Mreoow = +3.04 and Mrgiaw = +2.40, and an in- 
trinsic color of (mreo6w — Mrsiaw)o = 0.64. The BaSTI 
models—which have been updated recently, using new 
filter throughput curves released in 2020 October by the 
WFC3 team (Calamida et al. 2021)—yield similar ab- 
solute magnitudes of +3.01 and +2.38, and an intrinsic 
color of 0.63. Taking the averages, and using the ob- 
served baseline color of mre6oew — MFsiaw = 2.365 (Ta- 
ble 1), we find a color excess of E(mreoew — Mrsiaw) = 
1.73. Based on the values of A)/Ay given for the WFC3 
filters at the PARSEC website (derived from the in- 
terstellar extinction curves of Cardelli et al. 1989 and 
O’Donnell 1994, with Ry = 3.1), this color excess cor- 
responds to an extinction of Apgja4w = 3.33 mag. Thus, 
using an apparent magnitude at baseline of mpgi4w = 


19.581 from Table 1, we find a true distance modu- 
lus of (m — M)o = 13.86, giving a linear distance of 
Dg = 5.941.3 kpc. The quoted error is the formal 
uncertainty (dominated by the uncertainty in the spec- 
troscopic log g), but there likely are additional system- 
atic errors, given the large amount of extinction, the 
assumption of a standard value of Ry, the high sensi- 
tivity to log g, and other sources. For example, if a low 
value of Ry ~ 2.5, as deduced by Nataf et al. (2013) for 
the Galactic bulge region, were used, the distance would 
increase by ~2.8 kpc. 

The nominal 5.9 kpc distance places the source on 
the near edge of the Galactic bulge. However, it cannot 
be ruled out that the source is actually located closer to 
the center of the Galactic bulge, at a distance of +8 kpc. 
The star’s high radial velocity is consistent with stars in 
the bulge, but not large enough to rule out disk mem- 
bership. 

Taking the spectroscopy, photometry, and derived dis- 
tance at face value, these parameters correspond to an 
angular diameter of the source of ~0.0038 mas. 





9.3. Finite-source Effects Are Negligible 
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Figure 19. Location of the MOA-11-191/OGLE-11-462 


source star in the logg versus effective-temperature dia- 
gram (red filled circle with error bars), using stellar param- 
eters from Bensby et al. (2013). Superposed are theoretical 
isochrones from two sources. The blue curves show PARSEC 
isochrones, calculated for the source’s measured metallicity 
of [M/H] = +0.26, for ages of 2 to 10 Gyr. The orange curves 
are BaSTI isochrones for [M/H] = +0.3 and the same range 
of ages. 


We can estimate the effect of the finite size of a small 
source on the light curve by expanding A(wu) around u = 
uo, giving 


A(u) ~ A(ug) + A’(uo) (uw — uo) + 
+5 4"(wo) (u— uo)”, (18) 


where the prime symbol denotes the derivative. If we 
average over the face of a spherically symmetric source 
star, the term with the first derivative cancels out, leav- 
ing the quadratic term as the first one that contributes 
to the difference between finite and point-like source. 
Averaging the latter over an extent 6 explicitly gives 


uot+o 
A" (uo) (u — uo)? du 


uo—d 


= : AY (tip) 8. (19) 


(AA)s(6) == 


For u < 1, one finds A(u) ~ 1/u and thereby 
A’ (u) ~ 2/u3. Consequently, we find the largest differ- 
ences between finite and point-like source for the small- 
est u. For a source of angular radius 6,, we find a source 


size parameter p, = 0,/0n, and if up > p,, we can ap- 
proximate the magnification difference by 
Ps 


A= [(AA)s (VoE= WP) an 


Px 


With the angular diameter of the source 26, ~ 
0.0038 mas and the angular Einstein radius 0¢ 
5.2 mas, we find p, = 0,/0— = 0.00037 while up = 
0.00271. Inserting into Equation (20) gives AA = 1.5 as 
compared to A(uo) = 369 at closest angular approach 
between lens and source, in close agreement with a 
numerical evaluation of the average of the exact mag- 
nification function A(u). With a blend ratio g ~ 20, 
comparing 2.5log[(A + g)/(1+ g)|] for A and A+ AA 
gives a difference of 4 mmag between finite source and 
point-like source, which is of the order of the systematic 
errors of the photometry. 

Accounting for the finite source size will thus result in 
a very small change in uo, which has very little effect on 
other parameters, evidenced by their practically identi- 
cal values for our uo,+ and up, solutions. Similar to 
the case of the long-duration microlensing event OGLE- 
2014-BLG-1186 discussed by Dominik et al. (2019), the 
measurement of the parallax parameter 7 comes from 
the wings of the photometric light curve, not from the 
peak. Additionally, none of the astrometric deflection 
measurements are close to the peak. With the finite 
size of the source making little difference, we can safely 
neglect any limb-darkening effects. 


10. NATURE OF THE LENS 


In this section we discuss the nature of the lensing 
object. of the MOA-11-191/OGLE-11-462 microlensing 
event. We consider the mass of the lens, discuss con- 
straints on its optical luminosity, and consider whether 
it is a single object or could be a binary system. 


10.1. Mass 


As shown by Equation (16), the mass of the MOA- 
11-191/OGLE-11-462 lens can be determined from the 
values of its angular Einstein radius, @p, and the relative 
lens-source parallax, 7p. In §8, we derived Op = 5.18 + 
0.51 mas and 7p = 0.089 +0.014. These yield a mass of 
Miens = 7.141.383 Mo. 

An object with a mass this large cannot be a single 
(or double) NS or white dwarf. It can only be a BH— 
or an ordinary star, or conceivably a binary (or higher 


I 


| 
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multiple) containing stars, BHs, and/or other compact 
companions. To distinguish between these possibilities, 
we consider the observational constraints on the lumi- 
nosity and binarity of the lens. 


10.2. Lens Distance 


The mass determination for the lens is independent of 
assumptions about the distances to the lens and source. 
However, to constrain the optical luminosity of the lens, 
we first do need to determine its distance. This can be 
obtained from Equation (2), using the value of mpgs = 
0.463 + 0.051 mas (from §8), and the distance to the 
source (from §9.2), Dg = 5.9+ 1.3 kpc. 

These values yield a lens distance of DL, = 1.58 + 
0.18 kpc. Note that the derived distance of the lens is 
only a weak function of the adopted source distance. If 
the source were at the distance of the Galactic center, 
the lens distance would only increase to ~1.7 kpc. 











10.3. Luminosity Constraints 
10.3.1. From the Baseline Magnitude 





Main-sequence stars with a mass in the range 7.1 + 
1.3 Mo have a range of J-band absolute magnitudes of 
My, = —1.5+ 0.3, based on the tabulation® of stellar 
properties compiled by Pecaut & Mamajek (2013). 

The lens distance of 1.58 + 0.18 kpc, from the previ- 
ous subsection, corresponds to a true distance modulus 
of (m — M)o = 11.040.2. Thus a main-sequence lens 
of the measured mass would have an unreddened ap- 
parent magnitude of I ~ 9.5. The J-band extinction of 
the source (from §9.2) is ~3.3 mag, which is likely an 
overestimate for the lens, since it is considerably nearer 
than the source and may suffer less extinction. Nev- 
ertheless, this amount of extinction gives an expected 
apparent magnitude of a main-sequence lens of about 
I ~ 12.8, or brighter if the extinction is lower. The star 
would be yet brighter if it has evolved off the main se- 
quence to a higher luminosity. Since the baseline F814W 
magnitude of the source (plus lens) is ~19.6 (Table 1), 
a main-sequence or brighter stellar lens is conclusively 
ruled out. 











10.3.2. Direct Limits from Final-epoch HST Imaging 


A much stronger constraint on its optical luminos- 
ity comes from the fact that the lens was not de- 
tected in our final-epoch HST frames. The parameters 
Oe = 5.18 + 0.51 mas and tg = 270.7 + 11.2 days im- 
ply that the proper motion of the lens with respect to 








6 Online version of 2021 March 2, at http://www.pas.rochester. 
edu/~emamajek/EEM_dwarf_UBVIJHK-colors_Teff.txt 


the source is 6.99 + 0.75masyr~!. At peak magnifi- 


cation in 2011 July, the separation between lens and 
source was ~0.01 mas. Therefore, at the epoch of our 
last HST observation in 2017 (6.1 years after the peak), 
the lens-source separation had increased to ~42.6 mas. 
At this large a separation, a luminous lens would cause 
the PSF of the source in our images to show clear signs 
of elongation. In order to search for such a distortion, 
we subtracted a “library” F814W PSF” from the final- 
epoch F814W images (in the individual un-resampled 
_flc frames). We saw no indication of elongation. 

In order to place the most stringent limit possible on 
the lens brightness, we fit the above-referenced PSF to 
stars of similar brightness and color to that of the source, 
and adjusted the shape of PSF to optimize the fit (again, 
in the un-resampled frame). We saw a ~1% variation of 
the PSF, which is typical of minor, “breathing” -related 
changes in telescope focus. 

Figure 20 shows stacks of the source in our 2017 data 
at successive stages of analysis and illustrates the pro- 
cess of searching for photometric evidence of the lens. In 
particular, panel (D) shows the subtraction with the im- 
proved PSF. For illustration purposes, we show in panel 
(F) the residual pattern expected for a lens that is 10% 
the brightness of the source. There is no such signature 
visible in panel (D). We quantify this as follows. 

Using the improved PSF, we modeled the source im- 
age as a superposition of two stars separated by 42 mas 
(about 1 detector pixel) along the implied direction 
of proper motion. We did this simultaneously in the 
4 x 250 s F814W images from 2017. We found a best-fit 
flux for the lens of —1% of the source flux. We per- 
formed the same fits to similar-brightness neighboring 
stars, and found a distribution between —2% and +1%. 
Finally, we added to each neighbor star a companion of 
2.5% at the presumed location, and were able to recover 
the added flux at a level of 2.5+1%. We would clearly 
detect a lens if it were there. 

We therefore conclude that the brightest the lens could 
be at the presumed location in 2017 is ~1% of the source 
flux. We also explored how bright a lens could be at an 
offset of 42 mas in any direction from the source. We 
were able to rule out any lens brighter than 2.5% of the 
source flux. 

A luminosity of ~1% of that of the source corresponds 
to an apparent F'814W or J magnitude of ~24.6. At the 
distance and extinction of the lens, the limit on its abso- 
lute magnitude is M; = 10.2 (corresponding to a main- 








7 From _https://www.stsci.edu/hst /instrumentation/wfc3/ 
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Figure 20. Panels showing 10x 10-pixel stacks in the reference frame of 4x250-s F814W images from the 2017 epoch, centered 
on the source. (A) The original f1¢ data with only the bright neighbor to the upper right subtracted. (B) Stack of the best-fit 
library PSF. (C) The original minus the library PSF. The subtraction is done in the -flc frame, then the pixels are stacked. 
(D) The original minus the optimized PSF. (E) The previous panel with a neighbor added to the previous frame that has 10% 
of the source brightness and offset by 42 mas. (F) The source+neighbor pixels fitted with a single-star PSF; the centroid moves 


up and to the right to account for the neighbor. 


sequence star of about 0.2 Me or less, or a white dwarf 
of ~0.6 Me with a cooling age older than ~10" years). 

These considerations leave a single BH—or a BH with 
a mass greater than ~3.55 Mo in a binary paired with 
another BH, or with a NS or faint star—as the only 
viable possibilities for the lens. 


10.4. Constraints on Binarity 


Having established that the MOA-11-191/OGLE-11- 
462 lens is an extremely faint, or non-luminous, object 
with a mass of 7.1 + 1.3 Mo, we now consider in more 
detail whether the lens could be a binary system rather 
than a single, isolated BH. 

MOA-11-191/OGLE-11-462 was a high-magnification 
event, making its light curve especially sensitive to 
lens binarity. As described above, the event was in- 
tensely monitored photometrically by several microlens- 
ing groups in search of planetary-companion signatures, 
and no manifestation of a binary lens was seen. 

In fact, lens binarity will strongly affect both the pho- 
tometric and the astrometric signature if the angular 
separation between the components is around the angu- 
lar Einstein radius, 0, and it therefore has to be either 
much smaller or much larger. In either case, such a bi- 
nary lens would produce a diamond-shaped caustic. 

For a close binary, this caustic would be located 
around the center of mass of the lens, and both the pho- 
tometric and astrometric microlensing signatures would 
be strongly altered if the source-lens trajectory passed 
over it or came close. Specifically, for a mass ratio g and 
a separation dO, the caustic extends to (Schneider & 
Weiss 1986; Erdl & Schneider 1993; Dominik 1999) 

cr 2q 2 

sf = (+o? d’. (21) 
Comparing this to uo ~ 0.00271 yields the constraint 
d < 0.07 for the equal-mass case, with less than 11% 














variation over the range q € [0.4,1]. With DLO, = 
8.2 AU, this implies a separation smaller than 0.6 AU. 
Considering a systematic noise floor of ~5 mmag, devi- 
ations from a point-lens light curve exceed this for the 
given uo ~ 0.00271 and blend ratio g ~ 20 if d > 0.022, 
so that our data imply an upper limit to the separation 
of ~0.18 AU. 

In such a case, since one of the components has to have 
a mass of at least 3.55+0.55 Mo, the merging timescale 
for such a system is S10” years (Carroll & Ostlie 2006). 
Given that there is no evidence of active star formation 
in this location, this seems unlikely. 

In the case of widely separated binary components, 
the observed properties of the lens refer to one of them, 
leading to the same values as by assuming that the lens 
is a single isolated body, so that the mass determination 
of the putative BH remains valid. 

With Og referring to the mass of the object identified, 
in the wide-binary case one finds the caustic extending 
to (e.g., Dominik 1999) 





sf = a (22) 
Considering that any companion should be dark as well 
as the noise floor of our photometric data, we exclude 
any companion more massive than 10 % of the identified 
lens within 230 AU. 

To further quantify the limits on a binary companion 
to the lens, we performed a Rhie et al. (2000)-type anal- 
ysis in which we simulate hypothetical binary-lens light 
curves with the same epochs and properties as the real 
data. We then fit those light curves with a point-lens 
model to determine the potential signal from a given 
lens companion. For this test, we use a subset of the 
data that provides good coverage of the overall light 
curve and that is sufficient to provide dense coverage 
of the peak, where a potential companion would have 
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the greatest effect (Griest & Safizadeh 1998). Specifi- 
cally, we use the OGLE and MOA survey data sets and 
the microFUN data from Auckland Observatory, Kumeu 
Observatory, Vintage Lane Observatory, the Wise Ob- 
servatory 0.46m, and the Weizmann Institute. Setting 
the requirement that the y? due to the planet is >50, 
we find that the sensitivity limits can be approximated 
with a piecewise function 


1.54 for logg > —0.4 
0.3413 log qg + 1.677 
|1og(diim)| = sie 
for — 4.8 < logg < —0.4 
0.04 for logg < —4.8. 
(23) 


Thus, for g = 0.1, this analysis also excludes companions 
with separations 0.4 AU <a, < 180AU. 


10.5. The Lens is an Isolated Black Hole 


To summarize this section, we have found that the 
lens that produced the MOA-11-191/OGLE-11-462 mi- 
crolensing event has a mass of 7.1+1.3 Mo. There are 
stringent limits on its optical luminosity, ruling out the 
possibility that the lens could be an ordinary star of 
that mass. Its mass is greater than possible for a white 
dwarf or NS, or for a binary pair of white dwarfs and/ 
or NSs. Even if the lens is a binary, at least one of its 
components must still be a BH—but only under unlikely 
circumstances could it actually be a binary at all. In the 
rest of our discussion, we will assume that the lens is an 
isolated, single BH. 





11. PROPERTIES OF THE BLACK HOLE 


Table 6 summarizes the inferred properties of the BH 
lens that produced the MOA-11-191/OGLE-11-462 pho- 
tometric and astrometric microlensing event. In this 
section we discuss some properties of the lens and its 
Galactic environment. 


11.1. Galactic Location 


The BH lies in the direction of, but is considerably 
closer than, the Galactic bulge. The fifth row in Ta- 
ble 6 gives the rectangular components of its Galactic 
position relative to the Sun. The BH is located in the 
Galactic disk at a distance from us of ~1580 pc toward 
the Galactic center, and about 45 pe below the Galac- 
tic plane. This places it approximately between the 
Scutum-Centaurus and Sagittarius-Carina spiral arms of 
the Milky Way, in the terminology of Reid et al. (2019). 
There are no conspicuous star-forming regions, young 
stellar objects, or SN remnants at this location. 


11.2. The Surrounding Stellar Population 





Figure 21. Color-magnitude diagram for a sample of stars 
in the Galactic disk at similar distances to that of the black 
hole, selected from Gaia EDR3 as described in the text. The 
population at this location contains main-sequence stars of 
a younger population, and older stars with redder colors. 


Most of the stars contained in the small FOV of our 
HST frames lie at considerably larger distances than the 
BH itself. In order to study a larger sample of the rela- 
tively nearby stellar population surrounding the spatial 
location of the BH, we selected stars from Gaia EDR3 
that satisfy the following criteria: they lie within 20’ of 
the BH on the sky, have trigonometric parallaxes with 
fractional errors of less than 10% that fall in the range 
0.55 to 0.70 mas, and are brighter than G = 17. 

Figure 21 shows the CMD of this population in the 
Gaia magnitude and color system. The CMD is afflicted 
by severe differential extinction, which blurs the stellar 
sequences. Nevertheless, it shows two prominent com- 
ponents: a moderately young and sparse population of 
main-sequence stars bluer than Gpp — Grp ~ 1.0, and 
a denser older component of redder stars. There is lit- 
tle or no evidence for an extremely young population of 
massive stars, at the BH’s location. However, with a 
transverse proper motion of ~ 45kms~! (Table 6), it 
would take only 0.2 Myr for the lens to traverse 20’, so 
it’s possible that the lens originated much farther away. 


11.3. Space Velocity, Population Constraints, Natal 
Kick 
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Table 6. Properties of the MOA-11-191/OGLE-11-462 Black Hole Lens 




















Property Value Sources & Notes* 
Mass, Mens TL4£1.3 Mo (1) 
Distance, Dr 1.58 + 0.18 kpc (2) 
Einstein ring radius, 6 5.18 + 0.51 mas (3) 
Proper motion, (fla, /ts) (—4.36 + 0.22, +3.06 + 0.66) mas yr! (4) 
Galactic position, (X, Y, Z) (—4, —1580, —45) pe (5) 
Space velocities, (V, W) (+3, +40) km s~* (6) 





“Sources & notes: (1) This paper, §10.1; (2) This paper, §10.2; (3) This paper, Table 4; 
(4) This paper, derived from model in Table 4; absolute proper-motion components in 
Gaia EDR3 J2000 frame; (5) Galactic position relative to Sun: X in direction of Galactic 
rotation, Y in direction away from Galactic center, Z perpendicular to Galactic plane 
toward north Galactic pole; (6) Space-velocity components relative to Sun, assuming zero 
radial velocity: V in direction of Galactic rotation, W toward north Galactic pole; U 
component toward Galactic center is undetermined due to unknown actual radial velocity. 


The sixth row in Table 6 gives the components of the 
BH’s space velocity relative to the Sun in the V and W 
directions (in the directions of Galactic rotation and to- 
ward the Galactic north pole, respectively). The third 
component of the space velocity, U, is undetermined 
since the radial velocity of the BH is unknown; changing 
the radial velocity to —100 or +100 kms~! changes V 
and W only by about +0.2 and +2.8 kms~, respec- 
tively. The relatively low value of V indicates that the 
BH is not a member of an extreme Galactic-halo popula- 
tion, and may be moving in a roughly circular Galactic 
orbit (unless the radial velocity is large). However, it 
does have a moderately high W velocity perpendicular 
to the Galactic plane, so the orbit has a modest inclina- 
tion to the Galactic plane. Apart from these consider- 
ations, there are few constraints on the age of the BH, 
except that it has likely had time to move far from its 
birthplace. 

Figure 22 plots the absolute proper motion of the BH 
(black filled circle with error bars, from the fourth row 
of Table 6) and compares it with the proper motions of 
the neighboring Gaia stars from the sample described 
above. The directions of Galactic longitude and lati- 
tude are indicated at the lower left, showing that the 
dispersions are higher in the longitude direction than 
perpendicular to it. Members of the “young” popula- 
tion from Figure 21 are plotted with blue points, and 
the older population with red points. The older pop- 
ulation, which has proper-motion dispersions of o,,1 = 
2.59 mas yr + and Cy,b = 1.59 mas yr +, is dynamically 
hotter than the younger population, which has disper- 


sions of 0, = 1.08masyr~' and o,,) = 1.05 mas yr7!. 








The BH itself is a prominent outlier, predominantly in 
b, relative to the older and younger populations (~4 and 
6c, respectively). This suggests it may have received a 
“natal kick” from the SN explosion associated with its 
birth—assuming that the BH arose from the surround- 
ing population. Relative to the mean proper motion 
of the younger population, the BH’s motion is offset 
by about 6.0masyr~', corresponding to a tangential- 
velocity offset of ~45 kms~!. However, considering 
that its radial velocity is unknown, we cannot exclude 
the possibility that the BH might simply be passing 
through the surrounding population and could conceiv- 
ably be significantly older. In a recent paper, Andrews 
& Kalogera (2022) have given an extensive discussion of 
the possible natal kick of MOA-11-191/OGLE-11-462. 


12. DISCUSSION AND FUTURE WORK 
12.1. The Formation of Isolated Black Holes 


As reviewed in §1.1, it is generally believed that only 
stars above ~20 Mo collapse to form BHs. To form an 
isolated ~7 Mo BH, such as the one discussed in this 
paper, either its progenitor was a single star that ex- 
perienced considerable mass loss through stellar winds, 
it was produced through fallback in a weak SN explo- 
sion (e.g., Fryer et al. 1999), or it was born in a close 
binary that lost mass through a common-envelope mass 
ejection. In the binary scenario, the binary evidently 
became unbound, presumably when its lower-mass com- 
panion itself became a SN, and its mass loss and SN kick 
detached the binary. 

The MOA-11-191/OGLE-11-462 BH’s inferred mass 
of ~7.1 Mo lies remarkably close to the peak of the 
Galactic X-ray-binary BH mass distribution (see §1.1). 
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Figure 22. Proper-motion components from Gaia EDR3 
for a sample of stars at similar distances to that of the black 
hole, selected as described in the text. Typical uncertainties 
are about +0.05 mas yr‘ in both coordinates. Blue points 
represent a younger main-sequence population with Gaia col- 
ors of Ggp—Grp < 1.0, and red points an older and dynami- 
cally hotter population. The directions of Galactic longitude 
and latitude are indicated at the lower left. The black point 
with error bars shows the proper motion of the black hole, 
which is a significant outlier compared to the surrounding 
population. 





This may suggest a similar evolutionary path with mass 
loss through close-binary interactions, but a larger sam- 
ple of mass measurements for isolated BHs is needed to 
distinguish these different formation scenarios. 


12.2. Future Work 


There is no known X-ray or radio source at the po- 
sition of MOA-11-191/OGLE-11-462. Considering the 
high extinction toward the BH, it is possible that it lies 
in a region of high density of interstellar matter. If so, 
the accretion rate of material from the interstellar mat- 
ter by the BH could conceivably be large enough for it to 
be detected in deep X-ray and/or radio observations, as 
discussed in §1.2, and these would be worth attempting. 

Finally, we note the potential of future facilities, such 
as the Very Large Telescope Interferometer (VLTT), 
Nancy Grace Roman Space Telescope, and the Rubin 
Observatory, in this area. The GRAVITY instrument 
with VLTI is currently being upgraded to further im- 
prove its sensitivity (Eisenhauer 2019), which will make 
it possible to measure 0g from spatially resolved mi- 
crolensed images (Dong et al. 2019) of a large number 
of microlensed sources in the Milky Way. Roman will 


detect and characterize several thousand microlensing 
events (Spergel et al. 2015; Penny et al. 2019); with 
the same aperture size as HST, its astrometric capa- 
bilities are expected to be similar (WFIRST Astrome- 
try Working Group et al. 2019). According to current 
plans, Roman will observe microlensing events over a 
~2 square-degree region in the Bulge with a ~15-minute 
cadence for several months each year; thus it will make 
thousands of photometric and astrometric observations 
of every long-duration event, leading to very high accu- 
racy in both kinds of measurements. With such data, 
Roman should detect numerous BHs, yielding exquisite 
determinations of their masses, distances, and velocities, 
and providing important insights to our understanding 
of the formation and evolution of BHs. Rubin’s wide- 
angle survey will provide deep, long baseline photome- 
try and consistent data reduction, making it ideal for 
the discovery and photometric characterization of mi- 
crolensing events (Sajadian & Poleski 2019), probing 
many different lines of sight (and hence stellar popu- 
lations) across the Milky Way and Magellanic Clouds 
(Street et al. 2018a). In the bulge, it will complement 
the Roman data by filling in gaps in the Roman survey 
cadence (Street et al. 2018b; Strader et al. 2018). Rubin 
discoveries may need high-resolution timeseries imaging 
follow-up for full characterization of the events. 

Note: After our paper was submitted, a study by Lam 
et al. (2022) that includes an independent investigation 
of MOA-11-191/OGLE-11-462 was posted on arXiv. We 
have not used any of the measurements or results from 
that paper in our analysis. 


Based in part on observations made with the NASA/ 
ESA Hubble Space Telescope, obtained at STScI, which 
is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5- 
26555. Support for this research was provided by NASA 
through grants from STScI. HST data used in this pa- 
per are available from the Mikulski Archive for Space 
Telescopes at STScI,° under proposal IDs 12322, 12670, 
12986, 13458, and 14783. 

This work has made use of data from the Euro- 
pean Space Agency (ESA) mission Gaia (https://www. 
cosmos.esa.int/gaia), processed by the Gaia Data Pro- 
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Gaia Multilateral Agreement. 


8 https://archive.stsci.edu/hst /search.php 
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